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The mathematical theory of differential games is used to study the 
structure of optimal allocation strategies for some time-sequential combat 
games with combat described by Lan chester-type equations of warfare. The 
work reported here primarily concerns the application of existing differential 
game theory (including some recent developments of the author on necessary 
conditions of optimality for problems with state variable inequality con- 
straints) for the determination of optimal time-sequential combat strategies 
in several Lanchester-type differential games of tactical interest. Two 
time-sequential problems of optimal fire-support allocation are considered: 
the first examines the determination of optimal air-war campaign strategies 
within the context of ground-war objectives, while the second examines the 
determination of optimal time-sequential fire distribution strategies for 
supporting weapon systems (here taken to be artillery or naval ship gunfire) . 

In previous work the author has studied optimal air-war campaign 
strategies in a generalization of a well-known tactical air-war differential 
game model. This work extended previous work by considering valuation of 
surviving air-war assets and temporal variations in both the effectiveness 
of counter-air operations and also the return from ground support. The 
differential game model discussed above has been criticized, however, 
because it does not evaluate air-war tactics within the context of ground- 
war objectives. Consequently, this model produces suboptimization. Hence, 
a model is considered which is considerably broader in scope, does consider 
land-war objectives, but (unfortunately) is considerably less susceptible 
to closed-form analytic solution. For this new model partial results 
concerning the determination of optimal campaign strategies are obtained 
and contrasted with those for the model previously discussed above. 
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The determination of optimal fire distribution strategies for 
supporting weapon systems is a major problem of military operations research. 



The author has studied some aspects of this problem in previous research. 
Motivated by a paper which recently appeared in the open literature, optimal 
fire-support strategies are studied in an attack scenario. The dependence 
of optimal fire-support strategies on combatant objectives is examined for 
a problem previously considered in the open literature: necessary and 

sufficient conditions on the functional form of the criterion functional 
(terminal payoff only) for the optimal fire-support strategies to be inde- 
pendent of force levels (at least for a certain range of values) are 
developed. This shows that certain results previously reported in the 
literature are, in general, not true. A variation of this problem is then 
considered. In contrast to the previous work reported in the literature, 
the attrition structure of the problem at hand leads to the optimal fire- 
support strategy of the attacker requiring him to sometimes split his 
artillery fire between enemy infantry and artillery (counterbattery fire). 
Numerical examples are given. 
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See list of references at end of main body of this report. 



1. Introduction. 



This report is the fourth in a series of reports which document the 
author’s research on the structure of optimal (time-sequential) allocation 
policies for tactical situations described by Lanchester-type equations of 
warfare. The report at hand contains Appendices C and D of the author’s 
previous report [10]. The results of these appendices were discussed in 
general terms in [10], but it had not been possible to publish the appendices 
themselves at that time. For the reader’s convenience, portions of the main 
body of the report [10] are reproduced here as they pertain to the 
work at hand. It should be understood that in Appendices C and D of the 
report at hand when reference is made to Appendices A, B, and E, it refers 
to work previously published in [10] . 

By considering several combat scenarios, further insights have been 
gained into such important questions as: 

(1) How should fire be distributed over targets? 

(2) How should targets be selected? 

(3) Do target priorities change over time? 

(4) Do force levels affect the optimal allocation strategies? 

(5) How does the number of target types affect the optimal 
allocation strategies? 

(6) Do conflict termination circumstances affect the optimal 
allocation strategies? 

(7) How are the optimal fire distribution/target selection 
strategies affected by the nature of the attrition 
processes? 

(8) What is the effect of logistics constraints on such 



policies? 



(9) How does the uncertainty and confusion of combat affect 
optimal allocation strategies? 

Specific objectives of this research have been the identification of decision 
parameters and the further development of general principles for target 
selection, regulation of firing rate, and allocation of military resources 
in dynamic combat situations. A further discussion of research objectives 
and their relationship to defense planning problems can be found in the 
author’s past reports* [8], [9], [10]. It should be pointed out that al- 
though the quantitative analysis of tactics may be considered to be a 
relatively new field of military operations research, there is wide-spread 
interest in such work (see [4], [6], [7], [8], [15]). 

Our approach has been to combine Lanches ter-type models of warfare 
with generalized control theory (both deterministic and stochastic optimal 
control, dynamic programming, and differential games) [9], [13] (see also 
[14]). This has been done by considering a sequence of concrete situations 
for which optimization problems were formulated and solved by application of 
the mathematical theory of differential games. The structure of optimal 
time-sequential fire allocation strategies has been studied by considering 
the solutions for specific problems and comparing and contrasting these. 

The organization of this report is to discuss results in general 
terms in the main body and to leave details for the appendices. 

2 . Guided Tour of the Appendices . 

In this section we summarize the work which is contained in the 
appendices and explain why this work was done. The results reported here 

* 

A comprehensive review of pertinent literature prior to 1973 in the fields 
of (a) Lanches ter-type models of warfare, (b) differential games, and (c) 
optimal control of Lanchester-type attrition processes is also to be found 
in 19]. 
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may be considered to be extensions of our earlier work [8], [9], [10], [13]. 
Moreover, the work at hand lays the foundations for more extensive work on 
the quantitative analysis of tactics and applications of generalized control 
theory to problems of military operations research. 

As mentioned above. Appendices A, B, and E have been published 
separately in [10]. The report at hand contains only Appendices C and D. 

In Appendices B and C we apply the theoretical results of Appendix 
A to the study of optimal air-war strategies. This has been over the years 
a problem of lasting interest for defense planners. In Appendix B we consider 
a generalization of the tactical air-war game of A. Mengel [2] (see also 
[1]) in which aircraft effectiveness may vary over time and a residual value 
for aircraft at the end of the campaign is considered. Optimal strategies 
are characterized for this general problem. A complete solution (i.e. , 
explicit determination of optimal closed-loop strategies) is given for cases 
of constant coefficients. These results are much more general than previous 
ones and provide much insight into trade-offs between various planning 
parameters. Preliminary results are also obtained for problems in which 
aircraft effectiveness in the air war is constant but the value of ground 
support changes over time. The latter is an attempt to reflect the influence 
of combat air support on the ground war and hence avoid suboptimization. Two 
cases of time -depen dent "returns 11 from ground support are considered: 

(a) linearly-decreasing-with-time returns from ground support, 

(b) exponentially-decreasing-^with-time returns from ground 
support. 

One criticism that has been made of the model considered in Appendix 
B is that it does not evaluate air-war strategies within the context of 
ground-war objectives. In other words, this model (unfortunately) produces 
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suboptimization. Thus, in Appendix C we present preliminary results for a 
model which considers the development of optimal air-war strategies within 
the context of land-war objectives. These preliminary results show that 
the outcome of the ground war may be a significant factor in the deter- 
mination of optimal air-war strategies. 

The research reported in Appendices B and C was undertaken because 
an important question for defense planners is what are appropriate missions 
over the course of a campaign for tactical air power. The answer to this 
question has far-reaching implications for Navy air forces (both carrier- 
based and land-based) and, of course, the Air Force. Recently, the USAF 
Studies and Analysis Group has been using quantitative methodology [15] in 
trying to answer such questions. By considering the results given in 
Appendices B and C one can begin to see the effects of the nature of the 
combat optimization problem on the structure of optimal air-war strategies. 

The determination of optimal fire distribution strategies for 
supporting weapon systems is a major problem of military operations research. 
This problem is particularly of interest to the military tactician so that 
he may have a clearer understanding of the circumstances , for example, under 
which enemy infantry should be engaged by a supporting weapon system (such 
as artillery) and those under which "counter-battery” fire is to be preferred. 
Furthermore, the study of such tactical time-sequential allocation problems 
is relevant to the Navy mission of fire support (both by ship gunfire and 
by carrier-based air) . 

In Appendix D we further consider the determination of optimal fire 
distribution strategies for supporting weapon systems. In this case it is 
appropriate to call the supporting weapon system artillery. We consider 
extensions of some recent work by Kawara [3] concerning optimal strategies 
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for supporting weapon systems in an attack scenario which is a variation of 
a model considered by Weiss [16]. We first examine for what class of 
criterion functionals (terminal payoff only) the optimal fire-support 
strategies are independent of the force levels and then examine the depend- 
ence of the optimal fire-support strategies upon the form of the combat 
attrition model by considering slightly different combat dynamics than 
those considered by Kawara [3]: we assume that the attacker’s artillery 

produces ’’linear-law 11 attrition (see [11], [14]) against both the defender’s 
artillery as well as infantry. The development of a complete solution to 
the latter problem has involved solution phenomena not previously encoun- 
tered in a Lanchester-type differential game: the dual variables are 

discontinuous across a manifold of discontinuity of both players’ strategies. 

The research reported in Appendix D was undertaken to obtain a better 
understanding of the dependence of the structure of optimal fire-support 
strategies upon model form. Previous work by Weiss [16] (see also [2]) 
and Kawara [3] had given the impression that an optimal fire-support strategy 
consists in always concentrating all fire on one enemy target type. From 
our previous work [11] on a one-sided Lanchester-type, time-sequential fire 
distribution problem in which each of two enemy target types undergoes a 
’’linear-law” attrition process, we knew that an optimal fire distribution 
policy could consist in splitting one’s fire between available target types. 
This was our motivation for the examination of other attrition structures 
in Kawara’s problem. Additionally, Kawara’ s recent work [3] stimulated our 
research on the influence of combatant objectives on optimal fire-support 
strategies because of Kawara’s remarkable result that (at least within a 
certain range of force levels) optimal fire-support strategies do not depend 
on force levels. We were able to show that Kawara had used the only type of 
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objective function that yields this result. Thus, in general, the result 
is not true. 



3. Summary of Research Findings (same as in [10]) . 

Here we summarize our research results. We have (at least partially) 
accomplished the following tasks which were suggested for future research 
in our previous report [9]: (a), (b2) , (e) , (g) , (il) , ( j 2) , and (k) . 

Results are organized under the following headings: 

(1) solution methodology for time-sequential combat games, 

(2) insights gained into the optimization of combat dynamics, 

(3) implications for defense planning. 

Items (2) and (3) differ in that the latter is a management -oriented digest 
of the practical implications of our research, whereas the former is oriented 
towards a technical audience. Further amplification of results and conclusions 
is to be found in the appendices. 

a. Solution Methodology for Time-Sequential Combat Games . 

Our research has produced the following results on solution methodology 
for time-sequential combat games. Specifically, we have accomplished the 
following: 

(1) developed theory of state variable inequality constraints 
for Lanchester-type differential games (This theory is 
essential for all problems with, for example, negativity 
restrictions on force levels), 

(2) demonstrated use of theory of state variable inequality 
constraints for solving time-sequential combat games by 
developing solutions to several specific combat optimiza- 
tion problems, 

(3) developed complete solution to a fire support differential 
game. This involved the following technical difficulties: 

(a) singular surfaces, 

(b) discontinuity in the dual variables, 
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(4) developed methodology for determining the functional form 
of the criterion functional (terminal payoff here) which 
yields optimal strategies that do not depend on force 
levels (at least for a certain range of force levels) , 

(5) concluded that computational methods must give consideration 
to structural properties of optimal strategies in idealized 
versions like those considered in this report. 

b . Insights Gained into the Optimization of Combat Dynamics . 

Based on our study of the optimization of combat dynamics using the 
mathematical theory of differential games, we have reached the following 
conclusions : 

(1) the structure of optimal strategies depends on the 
following factors: 

(a) criterion functional, 

(b) combat attrition model, 

(c) battle termination model; 

the dependence is complex, and future research should 
concentrate upon the examination of numerous structures 
of tactical interest, 

(2) optimal fire-support strategies depend on the nature of 
the attrition process produced by the supporting weapon 
system; when the supporting weapon system produces casu- 
alties at a rate proportional to the product of the numbers 
of firers and targets of a particular type (see Appendix D 
for mathematical formulation) , the optimal fire-support 
strategy of the supporting weapon system has the following 
characteristics : 

(a) it depends on the optimal strategy for the 
enemy* s supporting weapon system, 

(b) the optimal distribution of supporting 
fires may be 

(i) to concentrate on enemy 
infantry , 

(ii) to split fire to avoid 
overkilling , 

(iii) to concentrate on enemy 
artillery, 

(c) counterbattery fire is the optimal strategy 
during the early stages of an attack and 
destruction of enemy infantry is the optimal 
strategy during the final moments (unless the 
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ratio of defender infantry to defender 
artillery is "extremely" large, in which 
case the infantry is always engaged) , 

(d) a split of supporting fires between enemy 
infantry and artillery as an optimal 
strategy can only occur when enemy infantry 
have some effectiveness against friendly 
infantry, 

(3) optimal air-war strategies are different for the case in 

which ground-war objectives are considered than for that in 
which they are not, 

(A) optimal air-war strategies (see Appendix B for scenario) 
depend on the following factors : 

(a) residual value of surviving aircraft, 

(b) force levels (especially when a side is 
going to lose the fight for air supremacy: 
faced with a future loss of air power, one 
may abandon the fight for air supremacy 
and get what ground support one can from 
his planes before he loses them all) , 

(5) cases with variations over time in weapon system ef- 
fectiveness are best studied as extensions of constant 
coefficient cases . 

c. Implications for Defense Planning . 

In our research reported here we have studied idealizations of 
allocation structures that commonly occur in defense planning studies. 

After studying these idealizations in order to gain insight into the 
structure of optimal strategies in the complex real-world problem, we have 
reached the following conclusions as to considerations that should be 
brought to the attention of defense planners. These results should be 
kept in mind by practitioners who perform more detailed computer simulation 
studies. 

(1) The combat optimization problem should be thought of as 
consisting of three parts: 

(i) combatant objectives, 

(ii) conflict termination conditions , 

(iii) combat dynamics. 
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Optimal combat strategies depend on all three of the 
above. More basic scientific research should be done 
on all three, especially the first two. 

(2) Force levels do affect (either directly or indirectly) 
optimal combat strategies. 

(3) It may be quite dangerous to generalize optimal combat 
strategies developed for specific problems. At present, 
more research is needed on specific problems in order to 
develop an understanding of the qualifications that may 
be necessary to make about specific study results. 

(4) Optimal air-war strategies must be based on ground-war 
objectives. Suboptimization results when this is not 
done. This suboptimization is a serious problem, since 
it may lead to winning the air campaign but losing the 
war. 



4. Suggested Future Research Tasks . 

After performing the research documented in this report, we feel that 
the current state-of-the-art for applying the theory of differential games 
to time-sequential combat games is such that more significant results may 
be readily obtained in the future. Moreover, this previous research pro- 
vides valuable perspective for identifying what appear to be the most 
important research tasks to be considered next. In our opinion the most 
important task is to examine (and then understand) the influence of objec- 
tives on optimal strategies. Another important task is to study the 
structure of optimal fire-support strategies. 

Based on our past research experience we feel that there is much to 
be accomplished in the future. Specifically, we suggest the following as 
future research tasks: 



(a) Examination of the effects of campaign objectives (as 
reflected by the criterion functional) on the structure 
of optimal campaign strategies. In all our past research 
we have with one exception always considered a linear 
utility for survivors in the criterion functional. It is 
of interest to examine how the valuation of survivors 
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affects the structure of optimal strategies • Other 
criterion functionals that might be considered are 

( 1 ) nonlinear valuation of survivors 
(either the difference or ratio of 
functions which are, for example, 
concave, quasi-concave, etc.)* 

( 2 ) the value of losses rather than 
survivors . 

It would seem appropriate to begin such an investigation 
by considering the simplest problem possible. 

(b) Further general study of optimal fire-support strategies. 

The actual military operations must be analyzed and 
decisions identified. For various cases, the appropriate 
scenario and, consequently, Lanchester-type model would 
be developed, and the optimization problem studied. It 
is felt that at this time it is most important to study 
various specific problems in order to gain insight into 
the structure of optimal fire-support strategies. 

(c) Further study of optimal fire-support strategies for the 
model given by equations (64) of Appendix D. Some 
further analytic work remains to be done on the fire- 
support differential game model considered in Appendix 
D. It remains to consider the effect of annihilation 

of fire-support resources upon optimal time-sequential 
fire distribution strategies. Additionally, further 
computational work on this problem would be worthwhile. 

(d) Further study the determination of optimal air-war 
strategies within the context of land-war objectives. 

Based on preliminary results given in this report, it 
appears as though the outcome of the ground war is a 
significant factor* in the determination of optimal 
air-war strategies and that optimal strategies devel- 
oped for a model not considering land-war objectives 
need not be optimal for a model which does consider 
them. This proposed research would extend the results 
given in Appendix C of the report at hand. 

A more comprehensive (although somewhat dated now) discussion of 
suggested future research tasks can also be found on pages 65-67 of our pre- 
vious report 19] (see also pp. 10-11 of [10]). Work still remains to be done 
on the following tasks suggested there: (b) , (c) , (d) , (f) , (h) , and (j). 



This is not considered, for example, in TAC CONTENDER [ 15 ]. 
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APPENDIX C: Optimal Air-War Strategies Within the Context of Land-War 

Objectives 



1 . Introduction . 

As discussed in Appendix B, an important problem of military operations 
analysis is the determination of optimal air-war strategies. This problem is 
not only of intrinsic importance but also its solution has applications to 
force mix/force structuring analyses. It is a topic of current interest 
[8] and has been for over thirty years. 

In Appendix B we considered a differential game model for the determin- 
ation of optimal strategies for the employment of a supporting weapon system 
(here tactical air power) which we called the t actical air-war game . The 
constant coefficient version of this model may be stated as 



maximize minimize {v x(t f ) 
U V x 



v y y(t f ) + 



(Bux-Avy)dt) , 



with stopping rule: t^ - T = 0 , 



subject to: 
(air-battle dynamics) 



dx_ 

dt 



r 




(l-v)ay , 
(l-u)bx , 



( 1 ) 



and 



with initial conditions x(t=0) = x^ and y(t=0) = y^ 



x,y > 0 
0 < u,v < 1 



(State Variable Inequality Constraints), 
(Strategic Variable Inequality Constraints), 



where 



x(t) and y(t) are the numbers of X and Y aircraft, 
respectively, at time t , 



a and b are constant (Lanchester) attrition-rate coefficients, 



C-2 



r and s are replacement rates, 

v (v ) is the value of one surviving X(Y) aircraft, 

X Y 

A(B) is the (constant) value of one Y(X) aircraft flying 
ground support for a unit of time, 

u(v) is the fraction of total X(Y) aircraft which fly ground 
support missions. 



A further discussion of the formulation of this model is to be found in 
Appendix B. The model is shown diagrammatical ly in Figure 1. 




Figure 1. 



Diagram of Tactical Air-War Game. 
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We have discussed that a shortcoming of (1) is that it does not evaluate 
air-war tactics within the context of ground-war objectives. This is because 
the objective functional (or any minor variation of it) does not directly 
reflect the outcome of the land war. This model suboptimizes. Assumptions 
have been made which result in a very mathematically tractable optimization 
problem, but one which, unfortunately, does not consider the entire system 
and its objectives. We will now consider a (still idealized) model which is 
considerably broader in scope, does consider land-war objectives, but (unfor- 
tunately) is considerably less susceptible to closed-form analytic solution. 

We will only present a few preliminary results with further study left as 
future work. 

Let us consider a war between X and Y forces, each of which is 
composed of both ground forces (called infantry) and air forces. The ground 
forces are engaged in combat, while each side has to decide the optimal use 
of its air forces, i.e. how to optimally allocate these resources between ground- 
support missions and counter-air missions (see Appendix B for a further discus- 
sion). For the purposes of our analysis here, let us assume that each side 

k 

wants to maximize the net worth of survivors at the end of a war which ter- 

kk 

minates (for the sake of convenience) with the annihilation of one side s 
or the other’s ground forces. Mathematically, the problem may be stated as 
the following. 



Further considerations for the determination of decision criteria are discussed 
in [3] and [4 ] . 

More ’’realistic' 1 battle termination conditions along with associated analytic 
difficulties are discussed in [71 . 
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2 2 

maximize minimize { £ v x (T) - £ v y (T)} with T unspecified, 



U V k=l \ 

with stopping rule: 



k-1 V k 



(a) x 1 (T) = 0, 



or 



subject to : 
(battle dynamics) 



(b) y j ( T ) = 0, 

dx 

d tT = _a li y i " Va 12 y 2’ 

dx 2 

dT = r 2 ' (1-v)a 22 y 2 ’ 

dy l 

dt ’ b ll X l ub 12 X 2’ 
dy 2 

dT = S 2 - (1 ' u)b 22 X 2’ 



with initial conditions 



x. (t=0) = x? and y.(t=0) = y? for i = 1,2, 

l l i i 



( 2 ) 



and 



x^,x^,y^,y 2 ^0 (State Variable Inequality Constraints), 

0 £> u , v 1 (Strategic Variable Inequality Con tsraints ) , 

where 





is 


the 


number 


of 


X 


infantry 


(i.e. 


X^) at time t, 


x 2 (t) 


is 


the 


number 


of 


X 


aircraft 


(i .e . 


X^) at time t. 


similarly 


for 


y x (t) 


and 


y 2 (0 > 







a . . 
J-J 



is a constant (Lanchester) attrition-rate coefficient (the 
rate at which one Y unit can destroy X^) , 



similarly for b . . , 



u(v) is the fraction of total X^(Y^) forces (i.e. aricraft) which 
fly ground support missions, 

and other symbols are used in a similar fashion as those for probJem (1). 
This situation is shown diagraminat icall y in Figure 2. 




Figure 2. Model for Determination of Optimal Air-War Strategies 
Within Context of Ground-War Objectives. 

In the above model (2), the objective of each side is to "win" the 
ground war as "efficiently" as possible. We see that this model then evalu- 
ates air-war tactics within the context of ground-war objectives. We will 
develop some preliminary results for the determination of optimal strategies 
for (2). Unfortunately, it is a very hard problem, and its final resolution 



may require making some simplifying assumptions or using computational methods. 
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Of particular importance is to compare optimal strategies for (1) with those 
of (2) to determine the effects of suboptimization. It is of great interest 
to determine the adequacy of air-war strategies which have been developed 
in a model like (1) by evaluating them in a model like (2) . For reasons of 
such a comparison, we have included replacements of both sides 1 aircraft in 
(2). It should finally be pointed out that our multiplier conditions for 
differential games with state variable inequality constraints play a signifi- 
cant role in the developments of this appendix. 

2 . Development of Basic Necessary Conditions of Optimality . 



i * 1,2 and that r 2 ,s 2’ a ll ,a 12 ,a 22’ b ll ,b 12’ b 22 > Xt may be shown usin S 

the results of Friedman that the differential game (2) has value (see pp. 210- 

230 of [2]) and that a saddle point exists in pure strategies (see pp. 234-235 
of [2]). We now will develop necessary conditions of optimality for (2). 

Applying the results of Appendix A, we have that the Hamiltonian is 
given by 



It should be clear that we have in (2) above that v ,v ^ 0 for 

A , Y . 



i i 



= p i^ _a ii y i _va i2 y 2^ + P 2^ r 2~^ 1 ~ V ^ a 22 y 2^ 



+ q l ^ _b l l x l _ub 12 x 2^ + q 2^ S 2~ ( 1-u ) b 22 x 2^ ” y ( t ) x 2 + r '^ t ^ y 2’ O) 



where 





^ 0 for x^ = 0 , 



and 



= 0 for y 2 > 0, 



n(t) 



^ 0 for y^ = 0. 
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Because of the stopping rule we need not consider the state variable inequality 
constraints (SVIC's) x^,y^ ^ ® h ave adopted above the following corre- 

spondence between state and adjoint variables: 

state variable dual variable 



X . 
1 


p i 


for 


i = 1,2, 


y i 


q i 


for 


i = 1 ,2. 



It is also convenient to define 



2 2 

$ 1 (T,x(T),^(T),o 1 ) = l v x x.(T) - l v y yi (T) + o 1 x 1 (T), (4) 



i=l i 



i=l i 



and 



2 2 

* 2 (T,x(T), X# (T),o 2 ) = l v x x.(T) - l v y y.(T) + ^(T). (5) 



i = l i 



i = l i 



The a djoint system of differential e quations for the dual variables is 

1 P _I _ 9H_ , 

a Xl b n q i’ 



dt 



dp ? SH 

dt“ a"x~ = u * b 12 q l + ( 1- u i ’ c ) b 22 q 2 + 



dq ] 

dT 



an 

ay. a ll P l’ 



dq, 

dt 



Ifr = V * a l2 P l + (1 - v * )a 22 p 2 • n - 



( 6 ) 

(7) 

( 8 ) 

(9) 



The boundary conditions at t = T for the adjoint variables will be discussed 
below. 

When x 2’^2 > extrema -^ strategic-variable pair, denoted as 

(u*,v*), is' determined by the max-mln pri n ciple . Hence, we consider 



maximize minimize H ( t , ^,^,p=0 ,n = 0 ,u , v) , 

Osu^l O^v^l 



( 10 ) 



which by (3) yields 
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and 



maximize x 2 u ^ b i2 ^ -q l ^ _b 22 ^ _q 2^ ’ 
O^u^l 



minimize (-v) (a 12 P ] _ a 22 P 2 > • 
O^V^l 



Thus, for > 0, we have 



u*(t) = 



1 for b 12 (-q 1 (t)) > b 22 (-q 2 (t)). 



0 f ° r b i2 ( - _q l < ’ t ^ < h 22 ( '~ q 2 ( ' t ^ 



It may be shown that there are no u-singular subarcs in the solution of 



We also have for > 0 



v*(t) = 



1 for a 12Pl (t) > a 22 p 2 (t). 



0 for a^p^t) < a 22 p 2 (t) 



It may be shown that no v-singular subarcs exist. 

2.1. Necessary Cond itions of Optimality on Constrained Subarc for x . 

X 2 X 2 

On a constrained subarc on which x n (t) = 0 for t ^ t ^ t n 

2 dx 6 X 1 X. 

strategic variable v*(t) is determined by - = 0 for < t < 

and hence 

r 0 

(1-v*) = 



a 22 y 2 (t) ' 



or 



r 2 X 2 X 2 

v*(t) = 1 tT\ for t < t < t. . 

a 22 y 2 (t) e £ 

dx 2 X 2 X 2 

Since we must have 0 & v* ^ 1 and - — (t ) < 0 (here t denotes z 

dt e e 

hand limit), such a constrained subarc can only occur when 



X 2 X 2 
r 2 < a 22 y 2 (t) f ° r C e * c * 4 ’ 



( 11 ) 

( 12 ) 

(13) 
( 2 ). 

(14) 
the 

(15) 
lef t- 



( 16 ) 
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X 2 X 2 X 2 
because r < a y (t ). Since x_(t) = 0 for t & t £ t„ , the dimension 
2 22 2 e 2 e J6 

of the state space may be thought of as being reduced by one, and hence the 

3H 

dimension of the adjoint variable space is reduced by one through — = 0* 

O V 



This yields 



X 2 X 2 

a l2 P l^ = a 22 P 2^ f ° r C e s c * t#, • 



22 K 2 



(17) 



d 3H 

The multiplier y(t) is determined by — (— ) = 0 and hence 

dt dv 

p(t) = •^-Ja 22 {u*b 12 (-q 1 ) + (l-u*)b 22 (-q 2 )} - a 12 b 1;L 



( 18 ) 



X 2 X 2 

for t < t < t 0 . 
e Jc 

Thus, the necessary condition of optimality y(t) ^ 0 leads to that on a con- 
strained subarc with x^(t) = 0 for a finite interval of time we must have 

X X 

a 22* u * b l 2^ -q l^ + ^ 1-u *^ b 22^"" q 2^ ^ a 12 b ll^ _q l^ for fc e 2 < c < ‘ U9) 

X 2 X 2 

When y 2 > 0 and x (t) = 0 for t^ £ t t^ , we have that v*(t) is 

dx 2 

determined by the condition that =0. It should be noted that (17) 

implies that (12) is identically satisfied. In this case only u* is to be 
determined by the max-min principle (10). However, when x (t) = 0, then 
the max-min principle (or, equivalently, (11)) no longer yields that we must 
have u*(t) given by (13). However, we may still take (13) to hold by a con- 
tinuity argument in which we consider x^(t) = e > 0 and then let e 0 . 

Thus, (13) may be used to write (19) as 

X 2 X 2 

a 22 maximum ( b 12 ^" q i^» b 22^ _q 2^ fea i2 b ll^~ q l^ f ° r C e < c < * ( 2 °) 

It remains to discuss the corner conditions which must hold at entrances 



to and exits from constrained subarcs on which x (t) = 0 for a finite interval 

X 2 X 2 

of time. At entrance to constrained subarc with x^(t) = 0 for t ^ t ^ t 

z e x, 

(see Figure 3 below), we must have 
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X X X 

P i (t e } = P i (t e } = P i (t e > for 1 



= 1 , 2 , 



X X X 

= q i (t e > = q i (t e > f ° r 1 



= 1 , 2 , 



( 21 ) 



X X 

H(t e 2 ) = H(t e ), 



X XXX XXX 

where H(t denotes H(t 2 , x (t 2 ) ,x (t 2 ) =0, y,(t 2 ) ,£.(t 2 ) 2 )> 

6 6 X 6 <u 6 0 0 

X 2 X 2 

y=0,n=0,u*(t ),v*(t ) = 1). It is readily seen that 

e e 



X X X X X 

H <t e 2 ) = + P 2 (t e 2 )(r 2 -a 22 y 2 ( te 2 )} 



XXX 
- q l (t e 2,b ll x l (t a 2) + 



and that 



XXX XX 

H( te 2 ) = - <l-r 2 /( a22 y 2 (t e 2 »)a 12 y 2 (t e 2 )} 



XXX 
- q l (t e 2 )b ll x l (t e 2> + 



so that H(t Q 2 ) = H(t e 2 ) =P 2 (t e 2 )(t 2 -a 22 y 2 (t e 2 )) = Py < t/) < r 2 " a 2 2 y 2 < C e 2 3 ^ 
whence 



X X 

a 12 P l ^ C G > ’ a 22 P 2 (t a >' 



( 22 ) 



since we must have < a 22^2^e ^ ‘ We t ^ ius see t ^ iat th e corne r conditions 

(21) are automatically satisfied when (15) through (19) hold, and thus they 
yield no new information. 



on 




Figure 3. Entry to and Exit from Constrained Subarc. 

^2 ^2 

At exit from constrained subarc with x (t) = t ^ t ^ t^ , we must 



have 



Pi(V> 



X 2 X 2 

Pi (t £ > = Pi (t £ > for 1 



- 1 , 2 , 



X 2 X 2 X 2 
q i (t £ > = q i (t £ > ’ q i (t * > f ° r 1 



= 1 , 2 , 



(23) 



X 2 X 2 
H(t £ ) = H(t £ Z ). 



It is readily seen that 



X 2 

H(t i > 



P l (t £ 2 ){ - a ll y i (t * 2) 



1 - 



a 12 y 2^ C !l ^ 



a 22 y 2 (C t > 



XXX 
■'•l U t 2>b ll x l (t t 2) + q 2 (t l 2)s 2’ 



C-12 



and that 



X 2 

H <'l > 



X XXX 

P l (t £ 2){ " a H y l (t £ 2) - a l2 y 2 (t il 2)} + P 2 (t Jl 2)r 2 



x 2 x 2 x 2 

)b ll x l (t £ } + q 2 (t £ )s 2 ’ 



X, 



x n 



so that H(t^ ) 



«<t t > - 



X, 



a 12 P l (t t > 



a 22 P 2 (t * > 



(24) 



Again, we see that the corner conditions (23) are automatically satisfied when 
(15) through (19) hold, and thus they yield no new information. 



2.2. Necessary Conditions of Optimality on Constrained Subarc for y . 



Similar analysis yields for a constrained subarc on which y (t) = 0 

Y Y Z 

2 2 

for t £ t ^ t„ that we must have 
e l 



u*(t) - 1 - 



S 2 Y 2 Y 2 
b 00 x 0 (t) for c e < t'h * 



(25) 



22 2 



Y Y 

2 2 

with s 2 < t> 22 x 2 (t) for t £ ^ t ^ t^ . Additionally, we must have 



Y Y 

b 12 ( _q i ( c ) ) = b 22 (-q 2 (t)) for t g 2 £■ t ^ t £ . 



(26) 



and 



n(t) ' b^ [b 22 {v * a 12' > l +(1 - v * )a 22 p 2 , ‘ b 12 a n p l ! * 0 £or ^ " H* ' <27 



The latter condition may also be written as 



Y Y 

b 22 {v * a l2 P l +(1_V * )a 22 P 2 } ^ b 12 a ll P l f ° r C e " b " C £. * (28) 



or by use of (14) (which was developed from (12) for y^ > 0) 



Y Y 

2 2 

t> 22 maximum ( a 12 P l ,a 22 P 2^ k b 12 a ll P l f ° r C e < t < t n ’ ( 29 ) 
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The corner conditions yield that at entry to constrained subarc on 

Y„ Y„ 

2 2 

which y 0 (t) = 0 for t £ t ^ t 0 , we must have 
Z e x, 

Y~ Y + Y 

P i (t u 2 ) = P i (t Q 2 ) = P i (* e 2 ) for i = !> 2 > 



q i (t e 2) = q i (t o 2) = q i (t e 2) f ° r 1 = 1 » 2 » 



(30) 



b 12 ( ~ q l (t e 2)) = b 22 (_q 2 (t e 2)) * 



Also, the corne r conditions yield that at exit from constrained subarc on 

Y Y 

2 2 

which y£(t) = 0 for t^ £ t £ t^ , we must have 

— "h 

Y Y Y 

P i (t £ 2) = P i ( 4 2) = P i (t * 2) for 1 = L ’ 2 ’ 



Y Y Y 

q i (t * 2) = q i (t £ 2) = q i (t * 2) f ° r 1 = 1 » 2 » 



(31) 



b 12^ q d fc £. ^ b 22^ q 2^ t 2. ^ 



Thus, the corner conditions yield no useful information that we did not 
already know. 



2.3. Boundary Conditions for Adjoint Variables . 
There are two cases to be considered: 

( 1 ) 

( 2 ) 



x 2 (T) = 0, 



y^T) = o. 



Case (1): x^T) = 0 . 

In this case we use the function defined by (4) in the transver- 

sality conditions. It should also be clear that we must have y^(t=T) > 0. 
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The boundary conditions at t = T for the adjoint variables are then 



P x (t=T) 


■ \ + °i 


where 


°i 


is an 


unrestricted multiplier, 










i = 0 


for 


x ? (T) 


> 


o, 




P 9 (t=T) 


= v - V 

Y 1 


where 








z 






(33) 


z 






± 


' ^ 0 


for 


x 2 (T) 


= 


0, 




q^t-T) 


= -v 
















(34) 


1 






, = 0 


for 


y 9 (T) 


> 


0, 






q,(t=T) 


= " V Y + v 2 


where 


v 2 












(35) 


Z 


Y 2 2 




z 


‘ £ 0 


for 


y 2 (T) 


= 


0. 





We may also write (32) as 

p^(t=T) = a where a is an unrestricted multiplier. (36) 

3 * 

We determine a from the transversality condition H(T) + -^r~ =0, or 

H(T,x(T),^(T),^(T),^(T),y,n,u^(T),v*(T)) = 0. (37) 

Using (3), we then have 

(-a){a 11 y 1 (T) + v* (T)a 12 y 2 (T) } + (v x _v 1 > ( r 2 - (l~v* (T) )a 22 y 2 (T) } 



or 



ct = 



+ v y u*(T)b 12 x 2 (T) - (v y -v 2 ){s 2 -(l-u*(T))b 22 x 2 (T)} = 0, 

( v x -\ > 1 ){r 2 -(l-v*(T))a 22 y 2 (T)}+v Y u* (T)b 12 x 2 (T)~ (v y -v 2 > { s 2 _ (l -u * (T) >b 22 x 2 (T) } 

(a 11 y 1 ( T ) + v*(T)a 12 y 2 (T)} 



If the system is on a constrained subarc with x (t) = 0 for t ^ t ^ 

x 2 

with t g < T, then it is clear that we must have y 2 (T) >0 so that v 2 = 0. 
Thus, the boundary conditions (33) through (36) for the adjoint variables take 



( 32 ) 



( 

T 



the form 
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P 1 (t=T) = a 



P 2 (t=T) = v x - v 1 where £ 0, 

( t=T) = -v y , 

q„(t=T) = -v . 

2 



where a is an unrestricted multiplier, (39) 

(40) 

(41) 

(42) 



Since we have been on a constrained subarc with x 2 (t) = 0 for a finite 

interval of time, Y must be using the control (15) with < a 22^2^ t ^ 

X 2 

for t ^ t £ T, and, in particular, 



r 2 < a 22 y 2^' 



(43) 



We also must have by (17) 



a l2 P l < ’ t "’ T ) a 22 P 2^ T ^‘ 
Combining (39), (40), and (44), we obtain 



(44) 



a 12 a a 22 (v X 2 V l ) ’ 



where must satisfy k 0 and thus 



(45) 



V 1 <a 22 V X 2 - a 12 o) a0 ' 



(46) 



Recalling that = ^ and substituting (45) into (38), we obtain 



a = - 



V 2 



a l2 r 

ta n y i (T) + — 



(47) 



22 L 



a 22 y 2 (T)_r 2 



so that by (43) we have a < 0. It should be noted that (46) is automatically 
satisfied with > 0. We also have 



P 2 a-i) = - 



12 



22 



a < 0 . 



(48) 
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Thus, we have proved 



THEOREM 1: 


If the system is < 


on a constrained subarc with x 0 (t) = 0 




x 0 




x 0 




for t 2 £ t £ T 


with 


2 

t < T, then we must have 




e 




e 




r 






(a) 


y m > — > 0 


(so 


that q 2 (T) = -v ) , 




22 




2 


(b) 


a < 0 given by 


(47). 





8W 



Recalling that a = p^T) = ^ ^ 



where W = W(t , ,x^ ,y^ ,y ^ denotes 



Bellman’s optimal value function [1], we see that there is a negative marginal 
return to X for having more infantry. This is understandable, since if he 
did have more forces (but still lost the ground war), the war would be 

prolonged to the advantage of Y who has complete air superiority. 

Y 2 

If the system is on a constrained subarc with y 0 (t) = 0 for t ^ 

Y 1 e 

2 

t ^ T with t < T, then it is clear that we must have x^(T) >0 so 
e l 

that = 0. Thus, the boundary conditions (33) through (36) for the adjoint 
variables take the form 



P 1 (t=T) 


= a where a is an unrestricted multiplier, 


(49) 


P 2 (t=T) 


= V v 


(50) 




2 




q 1 (t=T) 


II 

1 

< 


(51) 




i 




q 2 (t=T) 


= -Vy + v 2 where \)^ £ 0. 


(52) 



Since we have been on a constrained subarc with y£(t) = 0 for a finite 
interval of time, X must be using the control (25), and, in particular, 
at t = T we have 



b 22 x 2 (T) ' 



u*(t=T) = 1 - 



(53) 



C-17 



We also must have by (26) 



b 12 ( ' q l (t=T)) = b 22 ( " q 2 (t=T)) 



(54) 



Combining (51), (52), and (54), we obtain then 



V 2 b 22 (b 22 V Y 2 " b 12 V Y 1 ) * °* 



(55) 



Recalling that y£(T) = 0 and substituting (53) into (38), we obtain 



+ 



a = 



■x 2 ‘ 2 ■ bf 2 v Y 1 tb 22 x 2 (T) - S2> 



(56) 



whe 



re we must have s ^ < ^22*2^ S ° that a > 0. Thus, we have proved 



THEOREM 2: If the system is on a constrained subarc with y (t) = 0 

Y Y 2 

2 2 

for t ^ t ^ T with t < T, then we must have 
e e 



(a) 


K (T) > b > o 


(so that p (T) = v ), 




22 


2 


(b) 


b 22 V Y 2 * b 12 V Y 1 > 




(c) 


a > 0 given by 


(56). 



Case (2): y^T) = 0 , 

This case is similar to the previous one only with the roles of X 
and Y interchanged. 



3 . Introduction of the Switching Functions S^(t) and S^(t) . 

It is convenient to consider the switching function S^(t) defined by 



s u (t) = b 12 ( - q l (t)) ■ b 22 ( ~ q 2 (t)) 



(57) 
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Then, we may re-write the extremal control law (13) for the strategic variable 
u(t) as 



for y 2 > 0: 



u*(t) = 



1 for S u (t) > 0, 



0 for S (t) < 0, 



(58) 



and it follows from (57) and the adjoint equations that 
for > 0: 



dS u | P l (b 22 a 12 b 12 a ll ) 
dt ” = \ 



for v* = 1, 



(59) 



b 22 a ]2 p 2- b 12 a ll p l ^ p l (b 22 a 12- b 12 a ll > £ot v **°- 



and 



dS 

-r— = 0 when y 0 = 0 for a finite interval of time, 
d t Z 



(60) 



Similarly, it is convenient to consider the switching function S (t) 



defined by 



£> v (t) = a 12 P 1 (t) - a 22 P 2 (t). 

Then, we may re-write the extremal control law (14) for the strategic 
variable v(t) as 
for > 0: 



/ 1 for S v (t) > 0, 

v*(t) - | 

^ 0 for S (t) < 0, 
v 

and it follows from (61) and the adjoint equations that 
for x^ > 0: 



dS 
v 

dt 



(- q i)( a 22 b i2 a i2 b n ) 



a 22 b 22 ( ' q 2 )_a 12 b ll (_q l ) * ("V (a 22 b 12" a 12 b ll ) 



for u* = 1, 
for u* = 0, 



(61) 



(62) 



(63) 



019 



and 



dS 

■j— = 0 when x = 0 for a finite interval of time. (64) 



4 . Synthesis of Extremal Strategic-Variable Pair . 

In this section we will partially synthesize the extremal strategies 
in one of four cases. We hope to study this important problem more thoroughly 



in the future. 


There are four cases to be considered: 


(i) 


b 22 a 


12 > b 12 a ll 




(1) 


a 22 b 12 > a l2 b ll 




(2) 


a 22 b 12- & a l2 b ll 


(id 


b 22 a 


12' * b l2 a ll 




(1) 


a 22 b 12 > a 12 b ll 




(2) 


a 22 b 12 * a 12 b ll- 



Comments : 

(A) The quantity ^22^12 ma ^ cons idered as th e rate of destruction 

by of Y^s kill capability against (see also [5], [7]). 

(B) The quantity ^i2 a il ma ^ cons idered as rate °f destruction 

by X^ of Y^’s kill capability against X^. 



Case (II): b ^ 2 a 1 2 > b 1 2 a 1 1 and a 2 ^ 2 > a ^ ^ x . 

Let us first make some observations before we give a few preliminary 
results. For y 2 > 0 and p^(t) >0 by (59) we have 



dS 

dt 



H 



P l (b 22 a l2 _b 12 a ll ) > ° f ° r v * = 1 » 

b 22 S 12 P 2 ' b 12 a U P l * P 1 (b 22 a 12- b 12 a ll ) > 0 ^ V * “ °’ 



so that 
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dS 

dt 



^•(t) > 0 for > 0 and p^(t) > 0. 



(65) 



In a similar fashion it may be shown that 

dS 



(t) > 0 for > 0 and (t) < 0. 



( 66 ) 



Since we will synthesize extremals to (2) by working backwards from the end 
at t = T, it is convenient to introduce the backwards time t defined by 



t = T - t. 
dS dS 

It readily follows that . Thus, we may re-state the above 

results (65) and (66) as follows: 



(67) 



Also, 



PROPOSITION 


1: 


If 


y 2 > o 


and 


p x (t) 


> o, 


then 


dS 

—Si < 

dt 


0. 


> 


PROPOSITION 


2: 


If 


x 2 > 0 


and 


q x (t) 


< o, 


then 


dS 

v < 
dT 


0. 



Also, let us make some observations about the satisfaction of necessary 
conditions of optimality for a constrained subarc on which x„(t) = 0 for 

x 2 x 2 x 2 x 2 

t £ t ^ t with t < t. . Our assumption that a b 0 > a „b implies 

that for q 1 (t) < 0 we have ^~^i^ a 22 h 12 > ^ _q l^ a 12 b ll ‘ Thus » if b l2^ -q l^ 

£>b 22^ _q 2^’ then it: follows that a 22 b 22^~ q ^ a 22 b 12^~ q l^ > a 12 b ll^~ q l^' 

By the above, it follows that the two conditions a 22 b 12 > a 12 b ll anc * ^(O 

< 0 imply that 



a 22 maximum ( b 12 (-q^ ,b 22 (-q 2 ) ) > a i2 b n^ -q l^* 



Thus , we have proven 
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PROPOSITION 3: If q^t) 


< 0 and 


a 22 b 12 > a 12 b ll’ then the 


necessary 


condition 


of optimality (20) for being 


on a constrained subarc with x~(t) = 0 for 






x„ X„ 


2 

t ^ t £ 


2 

t ft with 


2 2 

t < t„ is satisfied. 


e 


£ 


e £ 



The significance of the above proposition is that it guarantees that (local) 
necessary conditions of optimality are automatically satisfied whenever the 
system is steered to = 0 and maintained that way for a finite interval 

of time. Similarly, we can prove 



PROPOSITION 4: If p^t) 


> 0 and 


b 22 a 12 > b 12 a ll- then the 


necessary 


condition 


of optimality (29) for being 


on a constrained subarc with y 0 (t) = 0 for 




Y. 




t 2 * t £ 


2 . . 
t. with 


2 2 

t < t n is satisfied. 


e 


i 


e £ 



Considering the above results it seems appropriate to examine under 
what conditions p^(t) < 0 and q^(t) > 0* From an integration of the adjoint 
equations (6) and (8), we have that for 0 ^ t ^ T 



p^(t) = p 1 (t=T)cosh» / a^ 1 b 1;1 (T-t) 



and 
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hi 



q^t) = q 1 (t=T)cosh/a 11 b 11 (T-t) - J p 1 ( t^sinh/a^b^ (T-t). (69) 



q 1 (t=T)sinh/a 11 b 11 (T-t), (68) 



Thus, we can state the following self-evident propositions. 



PROPOSITION 5: If p (t=T) > 0 and q^t-T) < 0, then 

P, (t) > 0 

and for all t t [ 0 , T ] . 

qj(t) < 0 
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Also, 



PROPOSITION 6: If p.,(t=T) ^ 0 and q 1 (t=T) < 0 (or, 

equivalently, p^(t=T) > 0 and q^(t=T)£0), 
then 

P, (t) > 0 

and for all t £ [0,T). 

q x (t) < 0 



Let us now consider the case in which Y wins , i.e. x^(T) = 0. 

There are two further subcases to be considered: 

(a) y 2 (t=T) > 0, 

(b) y 2 (t=T) = 0. 

In Subcase (a): y 2 (T) > 0 , the boundary conditions (33) through (36) 

for the adjoint variables become 



P 1 (t=T) = a 


where a 


is an unrestricted multiplier. 


(70) 


P 2 (t=T) ' V X 2 - V 1 


where 


= 0 for 
l £ 0 for 


x 2 (T) > 
x 2 (T) = 


o o 


(71) 


q 1 (t=T) = -v y , 










(72) 


q (t=T) = -v . 

2 










(73) 


us further assume 


that a 12 v X, 


> a v 
, 22 X„ 


and that 


b v > b v 

12 Y- 22 Y* 


If 



x 2 (T) > 0, then the assumption that b^ 2 v y > ^22 V Y ’ (57), (72), and (73) 

yield that S u (t= 0) >0 so that by (58) we have u*(t=0) = 1. The trans- 
versality condition (37) then yields that 



P x (t=T) = 



''x 2 t V (1 - V * (I))a 22 y 2 (I)1 + V Y 1 b 12 X 2 (I) - V Y 2 S 2 

{ a H y l(T) + v*(T)a 12 y 2 (T)} 



(74) 



C-23 



If we consider the case in which a ~p (t=T) ^ a v where p (t=T) is 

1Z ZZ X^ X 

given by (74), then S v (x=0) ^ 0 with q^(t=T) < 0 implies that S v (x) < 0 
for 0 < t < x , since then q.. (t) < 0 for 0 ^ t < so that by Proposi- 

ds v 1 Q 

tion 2 < ^ for 0 ^ x < Tq < x^. Also, if p 1 (t=T) ^ 0, then by 

Proposition 6 we have p^(t) > 0 and q^Ct) < 0 for all t C [0,T). Then 

Proposition 2 yields that S (x) < 0 for 0 < x ^ T, since S (t= 0) ^ 0 
dS V V 

and < 0 for 0 £ x ^ T. Thus, if 0 ^ a^^p^Ct^T) ^ a 22 V X * then 

v*(t) = 0 for 0 ^ t ^ T. Furthermore, since p (t) > 0, we have S (x=0) 

dS 1 U 

> 0 with — — (t) < 0 for all x ^ 0. Hence, the equation S (x = x ) = 0 
dx u 1 

has a solution, and there is a switch in X's extremal strategy at x - x^ 
in backwards time. To summarize, we have shown that for 

x (t=T) = 0, x 2 (t=T) > 0, y x (t=T) > 0, y 2 (t=T) > 0, 

and 

0 £a 12 p l (t=T) Sa 22 V X 2 and b 12 v Y 1 =■ b 22 V Y 2 - 
in backwards time the extremal strategic-variable pair is given by 



and 



u*(t) = 1 
v*(x) = 0 

u*(t) = 0 

v*(t) = 0 



for 0 t ^ 



for x < x T. 



Unfortunately, much more work remains to be done. 



In Subcase (b): y ^ (T) = 0 , let us briefly consider the further subcase 

1 Y Y 

2 2 

in which y„(t) = 0 for t £ t £ T with t < T. Then by Theorem 2 we 
2 e e 

must have 



x 2 (T) > 




> 0 , 



(75) 
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b 22 V Y 2 ^ b 12 V Y i ’ 



(t=T) = a > 0, 



and 



S 2 Y 2 

u*(t) = 1 - r 7 T for 0 ^ x < x 

b 22 X 2 e 



(76) 

(77) 

(78) 



By (31) and (57), we have 



S (t=t ) = 0. 
u e 



(79) 



Now (77) and q^(t=T) = -Vy < 0 imply by Proposition 5 that p^(t) > 0 
and (t) < 0 for all t £ [0,T]. We then have by Proposition 1 that 



dS 

dx 



— (x) < 0 for all x € [0,T]. 



(80) 



The conditions (79) and (80) yield that 

S (x) < 0 for all x £ (x Y ,T], 
u e 

so that by (58) we have 

Y 2 

u*(x) = 0 for all x £ (x ,T]. 



(81) 



It should be noted that by Proposition 4 the necessary condition (29) is 
satisfied on the constrained subarc. Thus, we see that this subcase can only 
arise when (76) holds and then 



u*(t) 



for 0 ^ t ^ T - x , 

e 



S 2 Y 2 
1 - j- x for T - x < t £> T. 

b 22 x 2 C ) 



Again, we must terminate our discussion with this preliminary result. 
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5 . Discussion . 

In this appendix we have applied the theory of differential games 
(including the theory of state variable inequality constraints (SVIC's) 
developed in Appendix A) to the study of optimal air-war strategies in a 
model which considers ground-war objectives. Our reasons for doing this have 
been explained in the introduction to this appendix. The results presented 
in this appendix are preliminary, being based on an initial examination. We 
feel that this is an important and promising area and propose such further 
work to ONR as a future research task. The complete determination of optimal 
strategies for (2) will be an involved task and must follow the solution 
methodology given in [6]. 

It should be pointed out that the model (2) may be thought of as a 
generalization of the supporting weapon system game of H. K. Weiss [9] (see 
also [6]). The model (2) reduces to this simpler case when r^ = = 0, 

a ll = b ll = °’ a l2 = a 22 = b 12 = b 22 = 1 ’ \ = \ = and V X 2 = V Y 2 = °* 



It should be noted that this corresponds to Case (II) above. Thus, study of 
the model (2) with r^ = = 0 and a^ = b^ = 0 is itself of considerable 

intrinsic interest, since it generalizes the results of Weiss [9] (see also 
[6]). It should be pointed out that the theory of SVIC f s has never been used 
before to study such a problem. 

From our preliminary results (which have been a characterization of an 
optimal strategy but not the actual determination of one) , it appears as 
though the outcome of the ground war is a significant factor in the determina- 
tion of optimal air-war strategies and that optimal strategies developed for 
the model (1) need not be optimal for the model (2). It would be of great 
interest to numerically compare the payoff to, for example, X from using an 
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optimal strategy for the model (1) in (2) with that from using the optimal 
strategy for (2). The sensitivity of the payoff to X's strategy is of 
particular importance. Thus, we propose to ONR to further study the model 
(2) and to compare results with (1). 

Our initial investigation here has produced an insight of particular 
importance for force-level planning. By Theorem 1 we see that increasing 
the force level (or the X ^ force level) late in the war may actually 

decrease the payoff to X (even if he follows an optimal strategy) if the 
system still reaches the end state of x^(t=T) = 0 and X 2 (t=T) = 0, since 
the war is prolonged in a losing cause to the profit of Y. Such phenomena 
should be examined more thoroughly in the future. This may be the most 
significant insight developed from our research. Intuitively, as was possibly 
the case in the war in Vietnam, if a side commits resources that are not 
sufficient to produce a favorable outcome in a time-sequential combat game 
but only enough to prolong the occurrence of an unfavorable outcome, then 
there may be negative marginal return from the resources, and it is foolish 



to use them accordingly. 
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APPENDIX D: Some Differential Game Models for the Study of Optimal Fire- 

Support Strategies . 

1 • I ntroduction * 

The determination of optimal fire distribution strategies for supporting 
k 

weapon systems is a major problem of military operations research. This 
problem in one form or another is probably one of the most extensively studied 
problems in both the open literature and also classified sources, although 
this observation apparently has not been made previously. The problem is of 
interest to the military tactician so that he may have a clearer understanding 
of the circumstances under which enemy infantry should be engaged by a support- 
ing weapon system (such as artillery) and those under which "counter-battery " 
fire is to be preferred. Furthermore, the study of this tactical allocation 
problem is relevant to the Navy mission of fire support (both by ship gunfire 
and by carrier-based air) . Related optimization problems (in which the support- 
ing weapon system is called "air") have been considered in the previous 
appendices . 

As we have noted before, Berkovitz and Dresher [3] make the point with 

kk 

which we agree that operational gaming is not a helpful device for solving 
a game or getting significant information about the solution." One must dis- 
tinguish between finding out how people make decisions and how they should. 

Most analysts agree that operational gaming is not a useful tool for answering 
the latter question. As Weiss has also emphasized [32], a model of an idealized 
combat situation is particularly valuable when it leads to a clearer understand- 
ing of significant relationships which would tend to be obscured in a more 

See [32] for a brief discussion of the distinction between a "primary" weapon 
system (or infantry) and a "supporting" weapon system. 

See [21] or [29] for a discussion of terminology and background. 
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complex (and "realistic") model. It is in this spirit that most of the work 
reviewed below has been done and in which we consider several simplified 
models for gaining insight into optimal fire support strategies. 

It seems appropriate to briefly sketch some of the work previously done 
on such problems and how our work here is related to this. The problem of 
the appropriate mixture of tactical and strategic forces (another aspect of 
the optimal fire support strategy problem) was extensively debated by experts 
during World War II. Some analysis details may be found in the classical 
book by Morse and Kimball (see pp. 73-77 of [20]). The problem was studied 
at RAND in the late 1940's and early 1950's (see [12]) and elsewhere (see [1]). 

It would probably not be too far-fetched to claim that this problem stimulated 
early research on both dynamic programming (see [2]) and also differential 
games (see [12], [14]). Today the problem of optimal air-war strategies is 
being extensively studied by a number of organizations (see, for example, [19], 
[23], [30]). 

In Appendix B we have considered a differential game model closely 

related to A. Mengel's "War of Attrition and Attack" (see pp . 96-105 of Isaacs' 

book [15]). Discrete- time versions of this problem for the determination of 

optimal "air-war" strategies (see also [12]) have been considered by a number 

of workers as time-sequential combat games [3], [4], [11] (see also [5], [8]). 

Another related problem was considered by Weiss [32], who studied the optimal 

* 

selection of targets for engagement by a supporting weapon system. Recently, 
Kawara [16] has studied optimal strategies for supporting weapon systems in 
an attack scenario which is a variation of the model considered by Weiss [32]. 

See [28], however, for a justification of the optimality of strategies determined 
by Weiss [32]. A general solution algorithm is also presented in this paper [28]. 



D-3 



Other recent work has considered various conceptual and computational aspects 
of time-sequential combat games [22], [23], [24]. 

Since our work here may be considered to be an elaboration upon and 
extension of Kawara’s fire support differential game [16], it seems appropriate 
to review the major results of his work and to relate our work here to it. 

In [16] Kawara considers combat between two heterogeneous forces, each composed 
of infantry (the primary weapon system) and artillery (the supporting weapon 
system). The problem is to determine the optimal strategy for each side in 
distributing the fire of its supporting weapon system over enemy target types 
according to the criterion of the force ratio of infantry at the end of the 
prescribed duration attack scenario. Kawara concludes that each side f s optimal 
strategy for the distribution of its supporting weapon system’s fire is to 
always concentrate all fire on the enemy’s supporting weapon system (counter- 
battery fire) during the early stages of battle (if the total prescribed length 
of battle is long enough) and then later to switch to concentration of all fire 
on the enemy’s infantry. He states that this switching time ’’does not depend 
on the current strength of either side but only on the effectivenesses of 
both sides’ supporting units” (p . 951 of [16]). Moreover, an optimal strategy 
has the property of always requiring concentration of supporting fires on 
enemy infantry during the final stages of battle. 

It is indeed curious that Kawara, who develops an analytic result which 
actually shows that the optimal strategy of one of the players (the defender) 
may depend upon force levels (see the expression for T” and its plot in 
Figure 4 of [16]), concludes that the switching time at which one changes 
his optimal strategy from concentration of all one’s supporting fires on the 

Kawara does not determine the optimality of extremal strategies developed for 
his problem (see [28] and also [25]). 



enemy 1 s supporting units to concentration upon the enemy's primary ones 
(i.e. infantry) is independent of force levels but depends only on the 
effectivenesses of both sides' supporting weapon systems [16]. Kawara's 
results as presented in his Figure 4 clearly show that the optimal strategy 
of one of the players (for his data set, the defender) is to switch (in forward 
time) from concentrating all supporting fires upon the enemy's supporting 
weapons (i.e. counter-battery fire in the case of artillery) to concentrating 
all fires upon his primary system at an earlier (forwards) time in those cases 
in which the player is going to lose all units of his supporting system before 
a critical time (here t = 0 - T^). We have seen that such a result occurs 
in the solution to the tactical air-war game considered in Appendix B (see 
Figure 3 of Appendix B) . Such behavior was apparently first noted by Professor 
Rufus Isaacs (see pp. 101-104 of [15]), whose discussion is more thorough than 
Kawara's [16] but who presents essentially none of the mathematical details 
in [15]. 

Thus, Kawara concludes in [16] that optimal strategies do not depend 
upon force levels. As noted above, this is only true provided that the appro- 
priate side's (in Kawara's numerical example, the defender) supporting weapon 
system is not reduced to a zero force level before a critical time. Let us 

assume, therefore, that neither side's supporting weapon system can be reduced 

& 

to a zero force level. Then the optimal strategies are independent of force 
levels, and this property of the optimal strategies is readily seen to be a 
function of the criterion functional for a given set of combat attrition 
equations. We will show that Kawara considered in [16] the only type of objec- 
tive function (i.e. criterion functional) which yields optimal strategies 

Initial force levels and the known length of battle may be sufficient to 
guarantee this for a given set (or range of values of) Lanchester attrition- 
rate coefficients . 
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independent of force levels. Kawara's criterion functional is a special case 
of a more general one considered in this appendix. The adoption of a criterion 
functional of this general form will be shown to be a necessary (as well as 
sufficient) condition for the optimal strategies to be independent of force 
levels. In other words, Kawara essentially used the only objective function 
that yields the above results reported by him. For any other form of criterion 
functional, the (forwards) time at which a side changes from concentration of 

all fire on the enemy's supporting units to his primary ones in following an 

optimal strategy will depend upon the force levels. 

The importance to military analysis of the above property of the optimal 
strategies is obvious. If the statement [16] that optimal strategies do not 
depend on current force levels is taken literally, it means that one may use 
an open loop strategy, i.e. u = u(t;t^,x^), for optimal results. This means 
that it is optimal to use a pre-programmed attack strategy (not depending upon 
(intelligence estimates of) initial force levels or even the dynamic evolution 
of combat). This is, indeed, a dangerous conclusion to reach, since it may 
mean the difference, for example, between a successful and an unsuccessful 
amphibious landing. Moreover, this conclusion is not even true for the problem 
considered by Kawara, and a minor variation in Kawara's model leads to optimal 
strategies which depend directly upon the force levels. 

Thus, our purpose in this appendix is to clarify the dependence of the 

structure of the optimal strategies (at least for the case in which neither 

As noted above, this is provided that the appropriate side's (in Kawara's case, 
the defender) supporting weapon system is not reduced to a zero force level. 

As qualified in the previous footnote, this is only true provided that the 
appropriate side's (in Kawara's case, the defender) supporting weapon system 
is not reduced to a zero force level. We have assumed above that this condition 
is satisfied. 
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side ! s supporting units are all lost) upon model form (i.e. the combination 
of the type of attrition equations and the nature of the criterion functional). 

We then consider a slight variation in Kawara’s problem for which the structure 
of the optimal strategy of one of the combatants is significantly different 
than that in the original problem [16]: the optimal strategy of one combatant 

depends directly upon the enemy f s force levels and is no longer to always con- 
centrate all fire on either the enemy’s primary or supporting weapon system. 
Furthermore, we will show that an optimal strategy in which a side divides 
the fire of its supporting weapon system between the enemy’s primary (infantry) 
and supporting systems can only occur when the enemy’s infantry has some 
effectiveness (in the sense of a non-zero Lanchester attrition-rate coefficient) 
against his infantry. The optimal strategy of one side to sometimes split its 
fire is very similar to that which occurs in a one-sided (optimal control) 
problem previously considered by us [26] (see also [27]) for the optimal dis- 
tribution of fire by a homogeneous force in combat against homogeneous forces. 

In [26] the enemy consisted of two weapon system types, each of which undergoes 
attrition at a rate proportional to the product of the numbers of firers and 
targets (referred to, for convenience, as ’’linear-law" attrition). In fact, 
this previous work of ours [26] was the motivation for our examination here 
of other attrition structures in Kawara’s problem. 

2. On the Relationship Between the Objective Function and the Structure of 
the Optimal Fire-Support Strategies in Kawara’s Fire Support Differential 
Game . 

In the recent paper [16] Kawara considered the attack of heterogeneous 
X forces against the static defense of heterogeneous Y forces. Both the 
X and Y forces are composed of two types of units: primary units (or infantry) 
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and fire support units (or artillery). The X infantry (denoted as X^) 
launches an attack against the Y infantry (denoted as Y^) . We consider 
that phase of the attack which may be called the "approach to contact. 11 This 
is the time from the initiation of the advance of the X^ forces towards the 
Y^ defensive position until the X^ forces actually make contact with the 
enemy infantry. It is assumed that this time is fixed and known to both 
sides and that infantry fire has negligible effectiveness against the enemy* s 
infantry during this time. During this time the fire support units remain 
stationary, and each unit has the capability to deliver either "point-fire" 
counterbattery fire against enemy artillery or "area fire** against the enemy* s 
infantry. 

It is the objective of each side to attain the most favorable force 
ratio of infantry possible at the end of the **approach to contact** at which 
time the force separation between the opposing infantries is zero and artillery 
fires must be lifted from the enemy*s infantry in order not to also kill friendly 
forces. Thus, the decision problem facing each commander is to determine the 
'*best** distribution of artillery fire over time between enemy infantry and 
enemy artillery in order to maximize the quotient of friendly infantry (numerical) 
strength divided by enemy infantry strength at the end of the approach to con- 
tact. This situation is shown diagrammatically in Figure 1. The reader is 
referred to Kawara’s paper [16] for further details of the model's development. 

It should be pointed out that this model also applies to the case of an amphib- 
ious landing and the determination of the optimal allocation of Naval ship-gun 



fire support. 
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1-u 




Figure 1. Diagram of Kawara's Fire Support 
Differential Game . 



Mathematically, the problem may be stated as the following, 

x i (t.) 



maximize minimi 
U 



. . . \ 
v lze \y 1 ( tf ) f 9 



with stopping rule: 



subject to: 
(battle dynamics) 



t f - T = 0, 
dx i 

i r = ~ va i x i y 2 » 

dx 2 

di (1 - ,)a 2>'2’ 

dyi h 

li ub l y lV 

dy 2 

Tt (1 - u)b 2 ,t 2 1 



( 1 ) 
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with initial conditions 

x^(t=0) = x° and y^(t=0) = y° for i = 1,2, 

and 

X 1 ,X 2’^1’^2 ^ ^ (State Variable Inequality Constraints), 

0 u,v 1 (Strategic Variable Inequality Constraints), 

where 

x^(t) is the number of X infantry (i.e. X^) at time t, 
x^(t) is the number of X artillery (i.e. X^) at time t, 
similarly for y^(t) and ^2^^’ 

a^ is a constant (Lanchester) attrition-rate coefficient 

(reflecting the effectiveness of fire against X^) , 

similarly for b^, 

and u(v) is the fraction of X(Y) artillery fire directed at opposing 

infantry forces. 

It should be noted that for T < +°° it follows that we will always have 
x^(t),y^(t) > 0. Thus, the only state variable inequality constraints (SVIC’s) 
that must be considered are * 2*^2 ^ 

Kawara f s results and conclusions [16] have been discussed in Section 1 
above. Let us assume that neither side's artillery is reduced to a zero 
force level . Then we will show that Kawara in his paper [16] chose only the 
possible type of criterion functional which yields optimal fire support 
strategies that are independent of force levels. We will do this for two 
cases: (1) only infantry considered in the criterion functional and (2) all 



force types considered in the criterion functional. 
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Before we develop the above results, it is convenient to first establish 
some preliminary results regarding the characterization of an optimal fire 
support strategy (i.e. the basic necessary conditions of optimality). For 
X l’ X 2’ y ] ,y 2 > the Hamiltonian is given by [7] 

= -P 1 va 1 x 1 y 2 - p 2 (l-v)a 2 y 2 - < l 1 ub 1 y 1 x 2 " q 2 ( 1_u ) b 2 x 2 ’ < 2 ) 

where we have adopted the following correspondences between state and adjoint 
variables : 



state variable 


dual variable 








x . 

l 


p i 


for 


i = 1 


,2, 


y i 


q i 


for 


i = 1 


,2. 



The adjoint system of differential equations for the dual variables is 



dP l . 


3H 


v*a ] y 2 p 1 , 


dt 


" 3x i ’ 


dp 2 _ 


9H 


u*b 1 y 1 q 1 


dt 


" 9x 2 " 


dq l _ 


9H 


u*b 1 x 2 q 1 . 


dt 


II 

\ — 1 

: ^ 
(TO 

1 


dq 2 _ 


3H 


V * a ] X 1 P 1 


d t 


’ 9y 2 ’ 


’ y l’ y 2 " 


0 , the 


extremal 



2 h 2 : 



2 f 2 ' 



(3) 

(4) 

(5) 

( 6 ) 



as (u*,v*), is determined by the max-min principle . Hence, we consider 



maximize minimize H(t,x^£,g^,c^,u,v) . (7) 

O^v^l 



Thus, for y ,y > we ^ ave 














( 1 


for 


s u (t) 


> o. 




U*(t) = ' 










(8) 




0 


for 


s (t) 


< o. 
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where the u-switching function, S u (t), is defined by 



s u (t) = b 1 (-q 1 (t))y 1 (t) - b 2 (-q 2 (t)) 



(9) 



Similarly, for x^,x^ > 0, we have 

V*(t) = 



1 for S (t) > 0, 
v 



0 for S (t) < 0, 



v 



( 10 ) 



where the v-switching function, S v (t), is defined by 



VO = a 1 p 1 (t)x 1 (t) - a 2 p 2 (t). 



(ID 



For future purposes, it is convenient to note the following: 



and 



p^(t)x^(t) = constant. 



q (t)y (t) = constant, 



( 12 ) 



(13) 



Both these results are readily established by direct computation using the 
state and adjoint equations. Using the boundary conditions for the adjoint 
variables at t = t^, one can easily show that the constant in (12) is 

* 

positive, while that in (13) is negative. Since an extremal trajectory is 
synthesized by working backwards from the end at t^ = T, it is convenient 
to introduce the backwards time t defined by 

T = T - t. (14) 

Then, it is readily shown that the switching functions satisfy the following 
differential equations 

dS 

= (l-v*)b 2 S v (x) - b 2 a 1 P 1 (t)x 1 (t) , (15) 

By an extremal trajectory we mean a path on which the necessary conditions of 
optimality are satisfied everywhere in time. 
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dS 

= (l-u*)a 2 S u (x) - a 2 b 1 (-q 1 (t))y 1 (t) . (16) 

dS dS 

From (15) and (16), it should be noted that - 7 —^- < 0 and < 0 for all 

dT dx 

t ^ 0. Thus, we see that singular subarcs (see [6], [12], [25]) are impossible 
for the problem ( 1 ). 

2.1. Case in Which Only Infantry Considered in Criterion Functional . 

We will first examine the case in which only infantry is considered 
in the criterion functional. It is convenient to re-write the problem (1) 
as follows: 



maximize minimize G(x 1 (T) ,y, (T) ) , 
U V 



with stopping rule: - T = 0, 

d x i 

— = -va lXl y 2 , 



subject to: 
(battle dynamics) 



dx 2 

5T ‘ -< 1 -''>a 2 y 2 . 

dy ! 

Tt ub i y iV 

dy 2 

Tt a-u)b 2 x 2> 



(17) 

(18) 



(19) 



with initial conditions 

x i (t= 0 ) = x° and y i (= 0 ) = y° for i = 1 , 2 , 

and 

X l’ X 2 ,y l’ y 2 ^ ® (State Variable Inequality Constraints), 



0 ^ u,v ^ 1 (Strategic Variable Inequality Constraints) . 



It seems appropriate to discuss the mathematical properties of the 
criterion functional (terminal payoff) G(x^ (T) ,y^ (T) ) . We assume that the 
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first partial derivative of G(x^,y^) exists for each of its arguments. Let 
us recall that G(x^,y^) is t ^ le P a y°ff to X. Thus, it seems reasonable to 
further assume that 



3 G 

3x 1 (T) 



> 0 



and 



3 G 

3y 1 ( T ) 



< 0. 



( 20 ) 



We now state and prove Theorem 1. It should be noted that all the 



results developed in the previous section are independent of the criterion 
functional. Let denote the "backwards time" of the first switch in fire 

distribution as determined by the optimal fire distribution strategies. Then, 



we have 



THEOREM 1: Assume that x 



2 ’ y 2 ’ 



and T are such that 



x 2 (T) > 0 



of 



x 1 (T) 



and 

and 



y 2 (T) > o. 



Then t is independent 



y 1 (T) 



if and only if 



G(x 1 (T),y 1 (T)) = F 



{x 1 (T)}^ 



y x (T) J 



where F(u) is a differentiable, strictly increasing 
function of u for u ^ 0 and k > 0. 



PROOF: (a) Since x^(T),y^(T) > 0 and we assume that x^(T),y^(T) > 0, 

the boundary conditions at t = T for the adjoint variables are 



' > i (t ‘ T) ■ > °- «i (t - T) 

P 2 (t=T) = 0, q 2 (t=T) 

Thus, by (12) and (13) we have 

p^(t)x^(t) = constant = x^(T) 



3y x (T) 
= 0. 



3G 



3x ] (T) 



> 0, 



0 , 



( 21 ) 

( 22 ) 



(23) 
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3 Q 

q 1 (t)y 1 (t) = constant = y^T) < 0 , 



(24) 



so that the differential equations for the switching functions are 

dS 



dx (L v*)b 2 S v (x) b 2 a l x i (T) D Xl (T) ’ 



and 



dS 



1 

9G 



dT ’ U-u*)a 2 S u (T) + a 2 b iyi (T )- 8 (T) , 



(25) 



(26) 



with initial conditions 



9G 



S u (,.o, --b iyi (T) 8 T) 



> o. 



and 



3G 



S v (t-0) - a^CT) 8](i(T) 



> 0. 



(27) 



(28) 



dS 



It should be noted that, for example, (l-v*)S (x) 0 so that - 7 -^ < 0 and 

dS V 

- 7 —^ < 0 for all x ^ 0 . 
dx 



(b) By (27) and (28) and the extremal strategies ( 8 ) and (10) , we 



have 



u~tx; = x 

7 for x [ 0 , x ] . 

‘ v*(x) =1 1 



(29) 



Using (29), the differential equations (25) and (26) are readily integrated 
to yield 



s u (x) = -b 1 y 1 (T) 9y / T) ~ b 2 a i x i (T) — (T) x, 



S V (T) ■ Vi (T) KpiT + a 2 b ! y i (T) a^Tfi T ' 



(30) 

(31) 



It is also convenient to write (30) and (31) as 
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and 



S u (T) 



x ] _(T) 



S v (T) 



-Hi 



Now let 





/■“S 

H 

'w'' 
i — 1 

L_! 


9G 


3G ( 


3y x (T) 


3x 1 (T) ) D 1 


x 1 (T) 


9G 


3x ] _(T) 




x l (T) 


9G 


3G ,f 


3x x (T) 


3 Yl (T) '| a l 


-y x (T) 


9G 


3y 1 ( T ) 


= minimum 


, u V 

( V T 1> 


5 



- b^T > , 



a 2 V 



(32) 



(33) 



(34) 



where is the solution to the linear equation S^(t=t^) = 0. 



(c) By (32) , (33) , and (34) , it is clear that is independent of 

x^(T) and y^(T) if and only if 



V T > 

y l (T) 9y_(T) 



= k = constant > 0. 



(35) 



It is convenient to write (35) as 



X 1 (T) 3x 1 (T) + ky l (T) dy^T) °* (36) 

The solution to the above linear first order partial differential equation 

(P.D.E.) is readily obtained by the method of characteristics (see p. 247 of 

[9]). Thus, G(x^ (T) ,y (T) ) is constant for 

dx x (T) dy x (T) 

"xj (T) “ ky~(T) 9 
or 

(Xj (T)} 

y x (T) 



constant . 



(37) 
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From (37), we see that 



G(x 1 (T) ,y 1 (T)) 



F 



{x 1 (T)} k '| 



(38) 



where F(u) is a differentiable, strictly increasing function of u for 
u > 0 and k > 0. Thus, Theorem 1 is proved. Q . E . D . 



Let t 2 denote the "backwards time" of the second switch in fire 

distribution as determined by the optimal fire distribution strategies. With- 

u v 

out loss of generality we may assume that . Th en we have 



u* ( t) = 0 
v*(t) = 1 



for t € (t ,t 2 ] . 



(39) 



We now state and prove Theorem 2. 



THEOREM 2: Assume that x^y^* anc * T are suc ^ that x 2 ^ > ^ 

and y 2 (T) > 0. Then is independent of x^(T) 

and y 2 (T) if and only if is. 



u v 

PROOF: (a) Let us first observe that the assumption implies that 

S (t=t 1 ) = 0 and S (t=t 1 ) > 0. 
u 1 v 1 



(b) Using (39) and (33) with t = t^, we may integrate (26) to 



obtain 



S v ( t) = 



' a 2 b 2' 



a x (T) ■■ ■ dp { — t 2 _ t 2j # 

1 1 9x 1 (T) a 2 b 2 1 
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(c) The "backwards; time" of the second switch in fire distribution 

ft 

x~ is the positive root of S (t=tJ =0 so that 
2 v 2 



T = /r_2__ _ T 2 

2 / a„ b „ 1 ’ 

whence follows Theorem 2- Q.E.D . 



(40) 



2 2 



As an immediate corollary to Theorems 1 and 2 we have 



COROLLARY 1 : 



Assume that x°, y ° » and T are such that X£(T) > ® 
and > Then the optimal fire support strate- 

gies for (17) through (19) are independent of force 
levels if and on ly if 



G(x 1 (T),y 1 (T)) = F 



{x 1 (T)> 



y x (T) 



where F(u) is a differentiable, strictly increasing 
function of u for u ^ 0 and k > 0. 



dS dS 

It should be noted that the fact that — — ^ < 0 and — — < 0 for all 

dx dx 

x > 0 guarantees that there will be no further switches in fire distribution 



as determined by the optimal fire distribution strategies. 



2.2. Case in Which All Force Types Considered in Criterion Functional . 

In this case we will show that there is no "simple” criterion func- 
tional such that optimal fire support strategies are independent of force 
levels. We consider the following differential game: 

ft u v 

It should be noted that the assumption that x^ < x^ implies that X 2 > T l* 
This follows from the facts (readily obtained from (32) and (33)) that x^ = 
x^ 1 = b^/^b^k) and x^ = a^k/(a^b^). The condition x^ < x^ then yields 
X-^ < 1 / 2 w ^ ence x 2 > T i vxa (^0) • 
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maximize minimize P (x (T) ,x (T) ,y (T) ,y » (T) ) , (41) 
U V J z I z 

with stopping rule and battle dynamics given by (18) and (19), respectively. 
We assume that the first partial derivative of P (x^ ,x 2 ,y^ ,y ) with respect 
to each of its arguments exists and that 



3P 

3x x (T) 



> 0 , 



_ 3P 
3x 2 (T)' 



> 0, 



3P 

3y x (T)' 



< 0, 



and 



3P 

3y 2 ( T) 



< 0. 



(42) 



We now state and prove Theorem 3. Let denote the "backwards time" 

of the first switch in fire distribution as determined by the optimal fire 

u v 

distribution strategies. Furthermore, we have that x^ = minimum 

where Tj 1 is the solution to the Linear equation S^(t=t^) = 0 and similarly 



for Then, we have 



THEOREM 3: Assume that and T are such that x 2 (T) > 0 

and y 2 (T) > both x^ and x^ are independen t 

of x^T), x 2 (T), y^T), and y 2 (T) if and only if 



P(x 1 (T),x 2 (T),y 1 (T),y 2 (T)) = R(A,B) , 



where 



A = 



x 1 (T) 



[y i (T)] 



1/k 



)k d-K)^y 2 (T) 



b l - 3 1 

b7 k IT x 2 (t) 

= y 1 (T)e e 



and 



K = k/k with k > k > 0. 

R(A,B) has the following properties: 

(1) the first partial derivative with respect to 
each argument exists, 



for 


A,B £ 


0 






(a) 


it is 
fixed 


strictly 

B, 


increasing in A 


for 


(b) 


it is 
fixed 


strictly 

A. 


increasing in B 


for 
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PROOF: (a) Since x^(T), y^ (T) > 0 and we assume that X 2 (T),y 2 (T) > 0, 

the boundary conditions for the dual variables are 



P l (tT) 3x,(T) 



9 P 

> 0, qj(t=r) = ^Yf) < °’ 



(43) 



P 2 (t T) 8 x 2 (T) > °’ q 2 (t T) 3y 2 (T) < °* 



Then, by (12) and (13) we have 



3P 



P 1 (t)x 1 (t) = constant = x ] (T) ^ ^ 



3P 



q 1 (t) y 1 (t) = constant = y x (T) ^ 



> 0, 



< 0. 



(44) 



As before, the differential equations for the switching functions S^( t ) and 
S^(t) are (25) and (26) with initial conditions 



and 



riP riP 

S u (,-0) - -b^tt) + b, yy^jy , 

3P 3P 

S v (t= 0) = a 1 x 1 (T) 9x ^ (T) - a 2 3 x 2 (t) * 

dS 



(45) 



dS 



It should be recalled that we have shown above that < 0 and < 0 



dx 



dr 



for all t ^ 0< 



(b) Our interest is in determining the class of criterion functionals 
for which the optimal fire support strategies are force level independent. 

We are also considering the special case in which x°, y^> and T are such 
that x 9 (T),y 2 (T) > 0. From (a) we see that S^(t= 0) ^ 0 yields that 
u* (t) = 0 for 0 ^ t ^ T and similarly for v*. Thus, for the purposes of 
our investigation here the cases in which S^(x=0) ^ 0 and/or S^x^O) ^ 0 
are not germane, since they are cases in which independence of strategies from 
force levels is not a meaningful property to study. Therefore, let us consider 



the case in which 
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S (t=0) >0 and S (t= 0) > 0, 
u v 



( 46 ) 



where the above are given by (45) . 

(e) By (46) and the extrema.] strategies (8) and (10) , we again have 
(29) so that the differential equations (25) and (26) with initial conditions 
(45) yield 



d? 



S u (T) “ b l y l (T) 3y. ] (T)' + b 2 8y 2 (T) - 2 u l"l vw 3 X] (T) 



3P 3p 

r - b n a,x, (T) 5 ~ tt T t, (47) 



and 



3P 



S v (T) " a j X l (T) 3 X]L (T)’ a 2 9 x 2 (T) ' 3y (T) 



9P + a.b.y. (T) ~ T t. (48) 



It is convenient: to write the above as 



V T) ~ X l ( f) 3~x. (T) 



3P 



2 3p 



' y l (T) 3y 3 (T) + b x 3y 2 (T) 



X 1 (T) 3x l (T) 



- b 2 ajt 



(49) 



and 



V T) 



' y l (T) 3y 1 (T) f 



3P 

X ] (T) 3x (T) 



2 3P 



a i 9x 2 (T) 



3P 



Y 1 (T) 3y.(T) 



3 2 b l T 



(50) 



(cl) By (49) we see that is independent of x^ (T) , x 9 ('i), y^(T), 

:nd y (T) if and only if 



3P 2 3P 

~y, (T) x — + 



l v ' 3y (T) b 3y (T) 

= — = constant > 0. 

3P k 



( 51 ) 



x, (T) 



3 X] (T) 
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Writing the above as 
X-j (T) 



3P 



3P 



3x, (T) + ky i (T) 3y (T) 



- - k 



9P 



9y 2 (T) 



= 0, 



r" * ~ J i' 

we obtain (again) by the method of characteristics [8] that 



(52) 



P(x 1 (T),x 2 (T),y 1 (T),y 2 (T)) = F(a,b), 



(53) 



where 



and 



a = 



{x 1 (T)> 



(54) 



b = y 1 (T)e 



h y 9 (T) 

b 2 2 



(55) 



(e) Similarly, is independent of terminal force level values 

if and only if P also satisfies the following first order P.D.E. 



x x (T) 



3P - 3P 2 3P 

3 X j (T) + Ry i (r) 3 Yl (T) a x 3 x 2 (T)“ 



(56) 



where k > 0. Since (in order that both and be independent of 

force levels) P must satisfy both (52) and (56), the latter condition 
places the following requirement on any F which satisfies (53) 



n r, 3F , - 3F 2 3F n 
(k-k)u — + kb — - — — =0, 
3u 3b a^ 3c 



(57) 



3P 



where we have used facts like y (T) ( . 

i dy ( 1 ; 



3F , 3F 

a 3a b 3b 



(This latter 



condition is merely a consequnece of the chain rule of partial differentia- 
tion.) In obtaining (57) we have made use of the fact that F = F(a,b,c.) 



with c = x 2 (T) also satisfies (53). 
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(f) Solving (57) by the method of characteristics, we obtain 



P(Xj (T),x 2 (T), yi (T),y 2 (T)) = F(a,b,c) = R(A,B) , 
k c 

a a 2 

where A = — - — , B = be , and K = k/k. Using (54), (55) and the fact 

b K " X 

that c = x^(T), we readily obtain the expressions for A and B in the 
statement of Theorem 3. 



(g) It is obvious that A,B ^ 0 . The conditions that 



9P 



> 0 and 



(T) 

> 0 yield that — - > 0 and — >0. It should be noted that R(A,B) 



must be such that 



3P 

< 0 and 



ap 



3A 



9y 1 (T) 



3y 2 (T) 



< 0. This implies that we must have 



3y 2 (T) 



< 0. Thus, we must have k>k (i.e. K>]), and the proof of 



Theorem 3 is complete. Q .E.D . 



Comment : An equivalent form of Theorem 3 is obtained by reversing the order 

of solving equations (52) and (56) . 



As before, let denote the "backwards time" of the second switch in 

fire distribution as determined by the optimal fire distribution strategies. 

u v 

Without loss of generality we may assume that and then (39) holds. 

We now state and prove Theorem 4. 
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THEOREM: 4: Assume that x^, y^. and T are such that 

x^CT) > 0 and y 2 (T) > 0 and assume that 
Let be independent of x^ (T) , x^CT), y^(T), 

and y^CT). Then is independent of x^(T), X 2 (T) 

y^(T), and y 0 (T) if and only if 



where 



P(x 1 (T),x 2 (T),y 1 (T),y 2 (T)) = S(A,B), 



( Xl (i)) k - K b7^ <T) - 

A = — - e 



{y 1 (T)} 



K+l 



and 



b y 2 (T) b t x (T) 
B = y^e Z o 1 Z 



k,K > 0. 

S(A,B) has the following properties: 

(1) the first partial derivative with respect to 
each argument exists. 



for 


A,B ^ 


0 








(a) 


it is 


strictly 


increasing 


in A 


for 




fixed 


B, 








(b) 


it is 


strictly 


increasing 


in B 


for 




fixed 


A. 









u v 

PROOF: (a) Let us first observe that the assumption ^ implies that 

S (t-t, ) = 0 and S (t=t 1 ) > 0. Furthermore, we have that 
u 1 v J 



s v (T=T i ) = h y i (T) 



9P 

ay l (T ) 





r /rp\ 9P a 2 9P 1 






V T) aXj^CT) a ± 9x 2 (T) 




| a i 


/ms 9P 


' a 2 b l T l 


L y i (T) 9y (T) j 





(58) 
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(b) Using (39) and (58), we may integrate (26) to obtain 



S V (T) = S v (T=T l ) f a 2 b L y l (T) 



9P 



»y, ( t) (T V 



a 2 b 2' 



a l X l (T) 9 ^TtT (, ‘ T 1 )2 ' < 55 



(c) The "backwards time" of the second switch in fire distribution 

* 



x 2 is the positive root of S (x^x^) = 0 so that 



DP 



T 2^ Q+ / a„b 



1 - 



2 3x 2 (T) 



2 2 



a l X l (T) 3x t (T) 



+ Q 2 - Tj , 



(60) 



where 



3P 



Q = 



3y 2 (T) 



Vl (T) 



DP 



It should be noted that (60) reduces (as it should) to (40) when 

<lx (1; 

— i p _ „ 0 

ay 2 (x) 



(d) The assumption that x^ is independent of x^(T), x ? (T), y^(T), 



and y (T) yields that it is both necessary and sufficient that P be given 



by (53). Rewriting (60) as 



3P 



1 DP 



+ T 2) 


i 1 - V2 1 


[ b x 9y 2 (T) b^ 9 x 2 (T) j 


CM 


°2 b 2 a l b 2 1 


1 x i (t) 31 L ( 



(61) 



Now, will be independent of x^(T), x^(T), y^ (T) , and y^(T) if and 

only if the right hand side of ((>1) is a constant. Hence, 

u v 

As noted in the proof of Theorem 2, it may be shown that x^ < x^ implies that 
x^ > x ^ . The proof follows that given previously above. 
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x l (T) 



3P 

9x 1 (T) 



k 3P - h _2_ 3P _ 
3x 2 (T)' k b L 3y 2 (T) 



(62) 



where k is a constant (which may be determined by ( 6 ])). Since P must 
be given by (53) , F must also then satisfy 



, 3F - 3F 

ka P kb — 

3a 3b 



k 3F 
bfTj 3c 



= 0, 



8 P 3F 

where we have used farts like x, (T) 7 — = ka — , etc, 

1 dx^ ( I ) da 

(63) we have made use of the fact that F = F(a,b,c) with 
satisfies (53). Solving (63), we obtain 



(63) 



In obtaining 
c = x 2 ^) also 



P(x 1 (T),x 2 (T),y 1 (T)) = F(a,b,c) = R(A,B) , 

K ^i 1 2° _ 

where A = a/b , B = be , and K = k/k. The expression for A and B 

given in the statement of Theorem 4 are readily obtained via (54) and (55). 



3P 



> 0 



and 



(e) It is obvious that A,B ^ 0. The conditions that f . 

oX, (I) 

3P 3 R 3 R ^ 

"z > 0 yield that vr > 0 and > 0. It should be noted that 

3x 0 (T) y 3A dB 

3P 



R(A,B) must be such that 



3y x (T) 



< 0 and 



3P 



9y 2 (T) 



< 0. The latter condition 



implies that we must have 



3A 



ay 2 (T) 



< 0. Hence, we must have K > 0, and the 



proof of Theorem 4 is complete. Q.E. D. 



Considering Theorem 3 and Theorem 4, we may conclude that there is no 
"simple” criterion functional which explicitly contains x^ (T) , x^(T), y^ (T) , 
and y£(T) and yields optimal fire support strategies for (41) , (18) , and 
(19) which are independent of force levels. We state this as Corollary 2 • 
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COROLLARY 2: Assume that x°, y£, and T are such that x 2 (T) > 0 
and y ^(T) > 0. Then if the criterion functional (Al) 
depends explicitly on x^(T), x 2 (T), y^(T), and y 2 (T), 
there is no "simple’' form for this criterion functional 
such that the optimal fire support strategies (for the 
differential game defined by (Al), (18), and (19)) are 
independent of force levels. 



3 . An other Model for Optima l Fire-Support A l location . 

In this section we will consider a slightly different version of 

Kawara's model. We will see, moreover, that the structure of the optimal 
fire support strategy for the attacker has a fundamentally different form 

than that for the previous problem (1) (i.e. Kawara's original problem [16]): 
the attacker must sometimes split his fire between the defender's primary and 
supporting units in order to "avoid overkill." Furthermore, the nature of 
this split depends upon the allocation of the defender's supporting fires. 

Let us again consider the attack of heterogeneous X forces against 
the static defense of heterogeneous Y forces. Each side is composed of 
primary units (or infantry) and fire support units (or artillery) . The X 
infantry (denoted as X^) launches an attack against the position held by 
the Y infantry (denoted as Y^) . Again, we will consider only the "approach 
to contact" phase of the battle. This phase is the time from the initiation 
of the advance of the X^ forces towards the Y^ defensive position until 
the X^ forces actually make contact with the enemy infantry. It is assumed 
that this time is fixed and known to botli sides. 

The X^ forces begin their advance against the Y^ forces from a 
distance and move towards the Y^ position using "cover and concealment." 
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The objective of the forces during the "approach to contact" is to 

close with the enemy position as rapidly as possible. Accordingly, small 
arms fire by the forces is held at a minimum or firing is done "on the 

move" to facilitate their rapid movement. It is not unreasonable, therefore, 
to assume that the effectiveness of X^ ’ s fire "on the move" is negligible 
against Y^ . We assume, however, that the defensive fire causes attri- 

tion to the advancing X^ forces at a rate proportional to the product of 
the numbers of firers and targets. Two possible justifications of this are 
as follows: because of the movement (and intermittent concealment) of the 

Xj forces and the distance involved, the Y^ defenders either (1) fire into 
a constant (but moving) area without precise knowledge of the consequences 
of their fire or (2) when they do aim fire at X^ targets, the time to acquire 
such a target is inversely proportional to the density of X^ forces and much 
greater than the time to kill an acquired target. Under each of these cir- 
cumstances it is well known that the postulated form of attrition occurs 
[6], [31]. 

During the "approach to contact," the fire support units remain 
stationary. Each unit has the capability to deliver counterbattery fire 
against enemy artillery or "area fire" against the enemy’s infantry. In 
other words, we assume that each side’s fire support units fire into the 
(constant) area containing the enemy’s infantry without feedback as to the 
destructiveness of this fire. On the other hand, the effectiveness of 
counterbattery fire is not symmetric with respect to the two combatants. We 
assume that the defender has the capability (for example, through the use of 
aerial observers) to sense when an enemy supporting unit has been destroyed 
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so that fire may be immediately shifted to a new target and that fire is 
uniformly distributed over the survivors. The attacker, however, either 
(1) does not have the capability to sense destruction of enemy fire support 
units accurately (and hence distributes his fire uniformly over the (constant) 
area occupied by the defender’s fire support units) or (2) if he does have 
adequate fire assessment capability, then target acquisition times (which are 
inversely proportional to the density of the enemy’s fire support units) are 
much larger than the time to destroy an acquired target. This leads to a 
attrition rate proportional to the product of the numbers of firers and 

Y 2 targets [6], [31]. 

It is the objective of each side to attain the most favorable infantry 
force ratio possible at the end of the "approach to contact" at which time the 
force separation between the opposing infantries is zero and artillery fires 
must be lifted from the enemy’s infantry’s position in order not to also kill 
friendly forces. Thus, the decision problem facing each side is to determine 
the "best" distribution of artillery fire between enemy infantry and artillery 
over time in order to maximize the infantry force ratio at the time of contact 
between the two infantry forces. This situation is shown diagramma tically in 
Figure 2. 



Alternatively, we may think that the attacker has massed so much artillery that 
X 2 targets are always easily acquired by Y 2 once an X 2 unit has been 
destroyed. Moreover, it will be assumed below that the initial X 2 force 
level is sufficiently large to guarantee that it is never driven to zero. 

This assumption is not essential for the structure of X 2 1 s optimal fire 
support strategy. A similar structural result may be obtained when X 2 f s 
attrition is the same form as that for Y 2 . We have made the above assumption, 
moreover, so that the resultant attrition model is most similar to Kawara’s [16] 
but yet yields significantly different results for the attacker’s fire support 
strategy. 



D- 29 







1-u 



Figure 2. Diagram of Fire Support Differential Game. 

The above assumptions lead to the following differential game with an 
attrition structure slightly different than that in Kawara’s problem [16]. 



maximize minimize 
U V 



{ 



VVi 

y x (t f ) ) ’ 



with stopping rule: 



subject to: 
(battle dynamics) 



t f - T = 0 , 

dt a ll X l Y l VJ 12 X i y 2’ 



dX 2 

— = -<l-v)a y , 



(64) 
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dy i 

di ub l y lV 



dy 2 

3T" * ■ (1 - u)b 2 y 2 x 2 > 



with initial conditions 

Xi(t=0) = x? and y^(t=0) = y° for i = 1,2, 

and 

X l ,X 2 ,y l ,y 2 (State Variable Inequality Constraints), 

0 ^ u,v ^ 1 (Strategic Variable Inequality Constraints), 

where all symbols are (essentially) the same as defined above for problem (1) . 

It should be noted that for T < 4<» it follows that we will always 
have x^(t), y^(t), and y£(t) > Thus, the only state variable inequality 

constraint (SVIC) that must be considered is x^ ^0. However, let us further 
as sume that the force level of the attacker's artillery is never reduc ed to 
zero . In other words, we consider the special case in which x° and T are 
such that X £(T) > 0* 



4 . Characterization of Optimal Fire Distribution Strategies for the Supporting 
Weapon Systems , 

It should be clear that in (64) above we have a il » ? > a 2 9^2 > 

Although the results of A. Friedman [10] concerning existence of value do not 

apply to the original version (64) of our fire-support differential game, they 

do apply to a suitably modified version. If we were to consider a version of 

dx 2 

this problem with = -(l-v)a^y 2 + r 9 where r^ > 0, then it may be shown 

(see pp. 210-230 of [10]) that this "modified" differential game has value and 



that a saddle point exists in pure strategies (see pp. 234-235 of [10]). 

We will now develop the basic necessary conditions of optimality for (64) . 
For x x’ X 2’^l’^2 > the Hamiltonian for (64) is given by 



H (t u > v > * "Pl (a n x i>'i +va J2 x i y 2 ) ” P 2 a 2 (1 " v)y 2 



- q 1 ub 1 y 1 x 2 - q 2 (l-u)b 2 y 2 x 2 , (65) 



where we have adopted the following correspondence between state and adjoint 
variables : 

state variable dual variable 



x . 
1 



y i 



p i 



q i 



for i = 1,2. 



The adjoint system of differential eq uations for the dual variables is 
d Pi 



dt 

dp 



1 

■ “llVl + v ‘ a 12 y 2 p l 



dt 

^1 

dt 

dq. 



3H 



3H 



= u*b 1 y 1 q 1 + (l-u*)b 0 y 0 q 



2 2 2 



3y, = a ll X l P l + u * b l X 2 q l 



3H 



with (T) = — , 



with P 2 (T) = 0, 



with q 2 (t> - - 7TT 



( 66 ) 

(67) 

, ( 68 ) 



(yp- 



dt 3y = v * a ]2 x l p l + 0~ v *) a 2 p 2 + (l _u *) b 2 X 2 q 2 With q 2^ = 



It should be recalled (see Appendix A) that H, g_(t) , and cj_(t) are con- 
tinuous functions of time except at manifolds of discontinuity of both U* and 
V* (see Section 4.3 below). 

When x^,x 2 ,y^,y 2 > 0, the extremal^ strategic-variable pair, denoted 
as (u*,v*), is determined by the max- min p rinc iple. Hence, we consider 



We use the term extremal to denote a trajectory on which the necessary condi- 
tions of optimality are satisfied at eve ry point in time. 



])— 3 2 





maximize minimize H(t 
O^u^l O^v^l 


>x.’)OR,’l> u > v > • 




so that 


u*(t) = ] 


j 1 for 
| 


S u (t) > 0, 


(70) 




1 


0 for 


S u (t) < 0, 




where the 


U-swi tching function , 


s u (t). 


is given by 






S u (t) = b 1 (-q 1 )y 1 




(71) 


and 


v*(t) = | 


r 1 for 


s v (t) > 0, 


(72) 






( 0 for 


S v (t) < 0, 




where the 


V-switching function, 


s (t), 

V 


i s given by 






s v (t) - 


= a il P l X l 


- a 2 P 2’ 


(73) 



It is readily shown that 



P 1 (t)x 1 (t) = constant = p^CDx^T) = —j , 



(74) 



ft' (q l y l } = a ll ~J y l (t) > °* 

y l 



(75) 



and 



dS 

dr = - a 2 ( i- u *)s u (t) - 



(76) 



We must further investigate the possibility of singuJar subarcs (see 
[26] or Chapter 8 of [7]). Let us first show that it is impossible to have 
a V-singula r subarc J In other words, v*(t) must be 0 or 1 almost 

We use capital letters to denote (closed-loop) strategies (e.g. U = U(t,x^,j^)) 
and corresponding lower case letters to denote the outcomes of strategies, 
or strategic variables (e.g, u(t) = U(t,x^,£)). 
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everywhere in time. The impossibility of a V-singular subarc is established 
dS 

by showing that > 0 for all t £ [0,T]. It is clear that 



(l-u*)S (t) ^ 0 for all t £ [0,T] . 



(77) 



Considering (75) and the fact that q^(T)y^(T) = - — <0, we see that 

y 1 

qfCtJyfCt) < 0 for all t t [0,T], whence follows the assertion via (76). 

3H 

It is possible, however, to have a U-singular subarc on which = 0 

(or, equivalently, S^(t) = 0) for a finite interval of time. There are two 
cases to be considered: (1) v* = 1 and (2) v* = 0. 



4.1. U-Singular Subarc on Which V* = 1 . 

When v* = 1, it is readily computed that 



and 



dS 

i 

dt 



x n 



V 



(a ll b l y r a 12 b 2 y 2 ) * 



(78) 



d 2 S 

i 

dt 2 " 



rx. 



X 2^ a n b i y l' >U * b l " ^ a i2 b 2 y 2^ 1-U *) b 2^ 



(79) 



3H 

Considering (71), the requirement that — = 0 yields the first condition for 

dU 



a U-singular subarc with V* = 1 



b i q l y l = b 2 q 2 y 2 ’ (8 °^ 

Considering (78) and (80) , the requirement that = 0 on a singular 

subarc on which — = 0 for a finite interval of time yields the second condi - 
du 

tion for a U-singular subarc with V* = 1 



a ll b l y l a 12 b 2 y 2 



(81) 
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On a subarc on which the first and second conditions for a singular subarc 

hold we additionally require that - 7 — 7 - f— 1 =0 so that (79) yields the 

dt v du 

singular strategic variable value required to keep the system on the singular 
subarc 



u*(t) = 



V b 2 



(82) 



^ ^ d 2 ^ 

Checking the generalized Legendre-Clebsch condition [17], [18] — f-i — 1 ^ 0, 

ou at v du y 



we find that on a subarc on which (80) and (81) hold we have 



fX, 



f 

_3r d 2 r 3Hi 
3u> 



^©1 ‘ “f ^ 2 ) 2 (a n (b j ) 2 y 1 + a 12 (b 2 ) 2 y ) > 0 . 

V 

We may write the equation of the U-singular surface (see p. 683 of [26] as 

(83) 



h _ a 12 b 2 

y 2 ' a ll b l 



for v* = 1 . 



4.2. U-Singular Subarc on Which V* = 0 . 

When v* = 0, it is readily computed that 

f 



dS u x x 

' ■ ~t (a ll b l y r a 12 b 2 y 2 ) - b 2 y 2 S v (t) - 
y l 



dt 



(84) 



and 



d 2 S dS 

-^T- = -u*b 1 x 2 + a 2 b 2 y 2 {-u*S u (t)+b 2 q 2 y 2 +p 2 x 2 [u*b 1 -(l-u*)b 2 ]}, (85) 



so that the first and second conditions for a U-singular subarc with 
are, respectively, (81) and 



a nVi 



a 12 b 2 y 2 + b 2 y 2 



x n 



{-s v (t)>, 



V* = 0 



( 86 ) 



This is a necessary condition for optimality. R. Isaacs [15] gives an 
equivalent condition (see [13]). 
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It should be noted (see [13]) that the above singular surface exists in 



^ space. It is convenient to write 

f 



1 



a 12 b 2 



Vl 



t(-S ( t ) } for v* = 0. 

f v 



y 2 a ll b l a 11 b ] x 1 
The singular strategic variable value is given by 



(87) 



u*(t) = 



b l+b 2 



x Vl 

P 2 V 



( 88 ) 



The requirement that u* ^ 1 yields that on a U-singular subarc with V* = 0 
we must have 

V” q 2^ y 2 & b l P 2 X 2’ 



It is readily checked that the generalized Legend re-Clebsch condition is 
satisfied . 



4.3. Discontinuity of Adjoin t Variables Across Manifold of Discontinuity of 
Both U* and V* , 

It is convenient to introduce the backwards time x defined by 

x = T - t. (90) 

From (82) and (88), we see that u*(r) must change, in general, discontin- 

uously from b 2^ b i +b 2^ to b 2^ b l +b 2^ 1_q 2 y 2^ P 2 X 2^ whenever v*(x) 

changes from 1 to 0. Let us consider the totality of trajectories on which 

this happens. The locus of points in the t,x^- space for such simultaneous 

switches is then a manifold of discontinuity of both U* and V* . Across such 

a manifold the adjoint variables need not be continuous (see Appendix A) . 

Let i = x (x,y) denote the backwards time at which v*(t) changes 
v v 

from 1 to 0. For future purposes, it will be convenient to consider a 
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simultaneous switch with u* changing from the singular control b^/ (b^-fb^) 
to 1. Then the manifold of discontinuity of both U* and V* is given by 



F(t,x,£) = t - T + t v (x,£> = 0, 



and 



= a n Vl - a !2 b 2 y 2 = °- 

Across the manifold of discontinuity, we have 



T, +. T, 3F 9G 

R. (%> " E. <%> - » ^ - 0 ^ 



(91) 



and 



T, +\ 

a. < T V > 



H(x + ) 

v 



T, 3F 3G 

a- <%> " p si " 0 »£ • 



3F . 3G 
H(T v ) + P 37 + ° 17 



or 



and 



3x 



T, +, T, 

P (t v ) = P (t v ) - 


V 

p 3x ’ 


(92) 


v» = ^l (T v }) + p 


9t 

3y x + 0a ll b l 5 


(93) 


V }) = (_q 2 (x v )) + p 


8x 

3y 2 ' ° a 12 b 2 ’ 


(94) 


H(xb = H(xj 


+ P • 


(95) 



Considering (71) and (73), it is readily shown that 



8t 



8t 



S U (T V } = a{a 11 (b 1 )2y 1 + a 12 (b 2 )2y 2 } + p(b l y l 3^7 - b 2 y 2 3^7 } * (%) 



and 



S v (T t> ’ - p(a 12 x l 



(97) 
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Recalling that u*(t ) = b /(b +b ), u*(t + ) = 1, v*(t ) = 1, and v*(t + ) = 0, 

v Z J. z V V V 

9l 

we may substitute (92) through (94) into (95) to obtain for a^xj, + 

3x 9 t 11 1 1 9x i 

a 2 y 2 + Vl x 2 3^ * 1 



P = 



a ll (b ] )Zy l X 2° 



3t 3t 3t 

(l_a ll X l y l TZ - a 2 y 2 ST - b l y l x 2 3^ 



(98) 



Then we may write 



S u (t v> ■ 0 { a U (b l )2 >’l + a 12 <b 2 )2y 2 



9t 9t 

a ll (b p 2y l X 2 (b l y ] 3^ ~ b 2 y 2 



(l- 



9t 



9t 



9t 



' iri J l 3 Xl a 2 y 2 3 x 2 b l y l x 2 3y^ 



I 



a n x iyi 



(99) 



and 



3t 



S (x + ) 

V V 



■ a U< b ] )2y l X 2 (a 12 X l 3 , 



- a 



9t 

A 

2 9x, 



-)o 



( 1 - 



3t 



3t 



3t 



a n x ,y 



n ri 3x, 



- a 0 y 



- b-y.x. 






2 J 2 9x 0 Tl 2 9y 



(100) 



5 . Synthesis of E xtrema l Strategic -Variable Pair . 

By the synthesis of the extremal strategic variable pair we mean the 
explicit determination (using the basic necessary conditions of optimality) 
of the time history of the extremal strategic variable pair (u*,v*) ' from 
initial to terminal time (see [15] and also [26]— [28]) . The basic idea is 
to trace extremals backwards from the terminal manifold (where boundary con- 
ditions for the adjoint variables are known) in such a way to guarantee the 

•f* 

It should be kept in mind that, for example, u*(t) = U*(t,x^ y) . 
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satisfaction of the initial conditions. Thus, it is convenient to introduce 
the backwards time t defined by (90) . 



5-1. Extremal Transitions in Strategic Variables . 

It seems appropriate to examine what are the possible transitions (or 

changes) in each strategic variable as we work backwards from the end (i.e. 

8S 

as t increases). It has been shown above that < 0 for all t £ [0,T]. 

Considering the boundary conditions (66) and (67) for the adjoint variables, 

it follows that S (t= 0) > 0. Thus 

v 



V*(t) = 



1 for 0 t x , 
v 



0 for t < t. 
v 



( 101 ) 



It will be convenient to refer to that phase of the planning horizon during 
which v*(t) =0 as V-Phase I (i.e. 0 ^ t < T - t^) and to that during 

which v*(t) = 1 as V-Phase II. 

Extremal transitions in u* for increasing t as shown in Figure 3. 
Thus, this figure shows what changes we might expect to observe in u* as 



we follow an extremal backwards from the end of the planning horizon at t = 0. 



During V-Phase II when v* = 1, 



dS 

i 

dx 



f f d 

b,x ,/y, > 0. When u* = 0, then — , 

111 , My _J 



dS 

i 

dx 



f 

'V 



V 

y f ' 

V 1 



^ a il b l y l" a l2 b 2 y 2' ) With S U (T=0) = 



<0. During V-Phase I when v* = 0 , 



y T 



^ a il b i y i” a l2 b 2 y 2^ + b 2 y 2 S v^ T ^’ Durin g both phases, the singular 



subarc may be exited with either u* =0 or u* = 1. Once u* becomes 0, 
it remains this way. The above statements will be further justified below. 
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Note: 



u = 4 
S 



b 2 1 




b l +b 2 J 


\ 


b 2 ] 




( b l+ b 2 j 





1 - 



q 2 y 2 



p 2 x 2 J 



for 



for 



v* = 0 , 



v* = 1. 



Figure 3. Extremal Transitions in u* for Increasing x. 
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5.2. Extremal Synthesis for x* < x . 
From the above we have 



so that by (70) we have 



S (x=0) 
u 




( 102 ) 



u*(x) = 1 for 0 x ^ x , (103) 

where x is the smallest zero of the equation S (x=x ) = 0. If the U- 
u u u 

singular subarc is reached in V-Phase II (see Section 4.1. above), then let 
us denote the backwards switching time at which u* changes from 1 to 
b 2 / ^^2^ aS T u # Clearly, it is necessary that x* < x^ for this singular 
subarc to appear in the solution. Thus, in general, there are two cases to 
be considered: 

(1) T U < v 

and (2) x* £ x . 

u v 

In the report at hand, we will only consider the former case with the latter 

following along the same general lines of development. Let us therefore, 

assume that a ^ , a ^ , , b ^ , b ^ , x^ , , y ^ , and y ^ are such that 

x* < x . We will give numerical results for this case below. 

—u v 0 



5.2.1. Extremals Near the Terminal Manifold. 



Recalling (78), we see that 



dS y 

“ — ~ > 0 (<0) if and only if — > (<) 

dx y 2 



a 12^2 

— . Considering (102), it is clear that S (x) > 0 for v* = 1 when 

a u b i u f 

y l f a l2 b 2 y i a 12 b 2 

> — . However, S (x) may change sign when — r < — - — — . The U- 

f a b 9 u J f a b 

y 2 11 1 y 2 11 1 



singular subarc occurs when both S (x=x*) = 0 and a^b.y, = a- 0 b 0 y 0 at 

u u 11 11 1 III 

x = x*. Thus, x* is the smallest root of 
u u 
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b l X 2 



■r + T* + 
r u 



. f f 

b l x 2 a ll>’l 



-b X X* 

1 2 u A 
e = 0 . 



(104) 



If y ± is given, then S u (t=t*) = 0 and = a 12 b 2 Y 2 may be combined 

to yield the value of required in order to reach the U-singular subarc 

f f f f b i (a nyf" b i x ?) 

(denoted as y *) . r fhus, for a y ± b.x. we have y„ = ~ — — . 

z 11 1 1 Z 2 b 0 /n , f .v 

2 a 12 (l-b 1 x 2 T*) 

(Other results are given below in Table I.) We denote the corresponding ratio 



, f f J 

of y^ to y as 






(y L J 



When S (t=t ) = 0 with a by < a 10 b 0 y 0 , it follows that x is 
uu II II 12 22 u 

the smallest root of the transcendental equation 

f 



b l - 



a H y i 



f i f 

. f , a li y i b l X 2 a u . 

' a i2 b 2 y 2 T u + ~T~' e " °* 



(105) 



It may be shown that 



9l u £ £ 

-r — > 0 where r = y,/y„. This latter result is useful 
3r J V 2 



in proving the following: 



THEOREM 5: 



T > T * • 


Then 


, u*(x) = 


V U 




f 


f 


v*(t) = 1 


for 


A 

r“H I *4— | 
>. 1 


[ y l] 

ly f i 

y 2 






y 2 



PROOF: The proof is by contradiction. Let r = yjVy^. 



for 


f (a) 


with 


>-2 



&)* 



In other words, we can find x such that S (t=t ) = 0 

u u u 



S (t) > 0 for 0 ^ t < x , 
u u 



( 106 ) 



for — > 



f y iV 



y 2 
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(b) Consider r = 



lx' 

f I f 

-2 
8t 



y 2 y 2 



+ e with e > 0 and such that t < t . 

u v 



3t 



Then it may be shown that > 0. In particular 



3r 



r= 



^ y i1* > This 



7r 



implies, however, that t > t* for r = r. 

u u 

u* = 1 



(c) Observe that \ . - for O^t^t <t so that y, /y„ = 

1 = 1 uv 12 



f, f b l X 2 T 
y l' y 2 6 * Hence > 



y l , » a l2 b 2 

^ ^ for r = r. 



(107) 



dS 



since then t (r=r) > t*. It has been shown above that - — > 0 for y n /y 0 > 
u N u dT J 2 

db 

a n 0 b 0 / (a n n b_ ) . Thus, (107) implies that ^ (t=t ) > 0, and hence 
1/ 2 II 1 dx u 



0 = S (t=t ) > S (t) for t £ (t -6 ,t ). 
u u u u u 



(108) 



This last statement (108) is a contradiction to (106) , and the theorem is 
proved. Q. E.D . 



Other results are obtained in a similar fashion. 



5.2.2. Field Construction . 

f f f f 

For a given set of terminal values x^ , , y^, and y^ an extremal 

may be traced backwards from the terminal manifold by a backwards integration 
of the state and adjoint equations combined with (70) and (72) (also (82) or 
(88)). By varying these terminal values the entire field of extremals (see 
p. 128 of [7]) may be obtained. 

The various types of extremals that may occur in the field of extremals 

are shown in Figure 4. This figure is representative of all our numerical 

results for t* < t (see Section 5.2.3 below). Pertinent information con- 
u v 

cerning each type of extremal is given in Table I. 



V-PHASE I \ / V-PHASE II 
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Figure 4. Identification of Various Types of Extremals for Which Information 
is Given in Table I. 
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Table I. Extremal Trajectories for Fire Support 



Problem with x < x . 

u v 



so 



yf'* 



( u * (t) = 1 * y i n 

< v*(t> - l u y‘ >y 2 f ' 



x is the smallest positive root of 
u 



1 * 
^ + x + 
f u 



b l x 2 



V f f 

b l x 2 a ll y l 



, I “b x T 

1 1 2 u A 

e =0, 



with the following bounds established: 



for a n y [ > b^: 



1 „ * , 1 
~f " T u " — f ’ 



a n y i 



b l X 2 



, f , f 

for a \\^ \ ~ b-^ x 2 : 



* 1 

x = 

u 



b l X 2 



p f . f 

for < b ! x 2 : 



1 „ * . 1 

, f T u f ' 

b l x 2 a ll y l 



for a 11 y 1 I* b^: 



f* _ b i (a n y r b i x 2 ) 

1 1 b f * 

2 a 12 (1 ’ b lVu ) 



f° r a lx Y^ = b i x 2 : 



f * a n b i y i 

'2 - a 12 b 2 



e. 



S v (r) 



= a^b . 



f , 


r f -i 


f M ( 


a l2 , a ll y l 


y i < 


a 2 b l (b^^x^) 2 . 



a n y i 

b l X 2 



X - 



a ll y l 
(t>l X 2 ) 2 



, f 

b iV 



Also 
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Table I. (cont.) - 1 



l. P 



II 

so : 



(concluded) 



Let 



* II 

S (t=t ). Also, on we have 

v u SO 



f , f 

f f b i x 9 x 

^(t) = x 1 exp(a 12 y 2 x + (e -1)} , 

b l X 2 



x 2 (t) = x 2 , 



, f 

f X o T 

/ \ t 1 2 
y i ( ^ ) = y 1 e 



y 2 (-0 = y 2 , 



and 



Pj (t) 



P 2 (t) 



qj_(T) 



q 2 (t) 



1 r r £ a il Y l , b l X 2 X 

-J exp {- [a^ 2 y 2 T + f~ “ 1 )]J 



b l X 2 



= b. 



f x ll 


/ 


fj 


{ 

\ 



a 



{- 



i - 



ii 



f 

a v ^ 

lr 1 



b l x 2 



T + 



f u f 

a ll y l , b l X 2 T 
2 — (e 

( b lX 2 ) 2 



b l x 2 



11 



b l x 2 



. f 
" b l X 2 X 



} 



= -a 



12 



f 






1)} , 
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Table I. (cont.) - 2 



2. P 11 : 

SI 



u*(t) = b 2 / ( b 1 +t> 2 ) 
v*(x) = 1 



for 

u v 



* II 

where t is determined in 1. On we have 

11 si 



SI 



S u (t) = 0, 



and 



a il b i y i = a i2 b 2 y 2* 

^ A 

is the smallest positive root of S v (t=t^) = 0, where 



ii* n* 

S v (T) = S V + a 2 b ! y l 



’ll 

y,' ' (0x b ) 2 



u*^ y 
q l , 

_ X 



III 

+ fl f 

6x ^J 



(T_T u } " 77^7 ^Plex^Cx-T*) ) } , 
(0x 2 ) 



with 0 = b^b^/ (b^+b 2 ) . An upper bound on t v is given by 



T V = a !2 /(a 2 b l ) - 



Also, on PgJ we have 



* 1 (t ) = xj exp | 



a ii y i +a l2 y 2 



ex. 






x 2 (t) = x 2 , 

f * 
0X (x~T ) 

yjd) ■ yf . 2 “ 

, , £ 8,, 2 < '-^ ) 
Y 2^ t) = y 2 6 



and 



u * ( 

Pi<t) = exp -J - 



u* f- 

a il y l +a i2 y 2 



ex. 



f ex,( T -T*) 

b u -i 



.]} 



u* 



1— 1 


fl) 


u* 




X 1 


V 




’ll 





f 


_ 


u* 


f y !^ 


. a ll 


q l 


f 


+ f 




lx l J 


0X 2 I 



with p 



, * a il 0X 2 (T " T u ) 

( T - T u ) ?— e 

(6x1)2 



}• 



f * 

-e x (t-t > 

2 u 



(T) “* - 0X 2 (T - T u ) 

q 2 ( T ) = q 2 e 






f * 
-6x 2 (t-t ) 



y i 0x 2 



y 7 0X 2 



12 /. - 0x 2 ( t-t *) 
“P " 6 



* 
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3. P 11 : 

A1 



Table I. (cont.) - 3 



/ u*(t) = 1 y f ryf|* 

< for 0 £ t £ t with -— > 1 Al 



v*( T ) = 1 



y 2 y 2 



t is the smallest positive root of S (t»t ) = 0, where S (t ) 
v ^ v v v v v 

is given in 1. An upper bound on is given by 

: a i2 



V a 2 b l ' 

It has been shown that S (t) > 0 for 0 £ x s. t . The solutions to the 

u v 

state and adjoint equations are the same as those for given above. 



r ii- { 



u*(t) = 1 

v*(t) = 0 



for T & T £ T 
V u 



We have that S (t) <0 for t > t and that 
v v 



df (T) = b 2 y 2 S v (T) + 



“£ (a ll b l y l (T) - a 12 b 2 y 2 ) - 



Also, on P^ we have 
dx. 



^r = a n x i y i with x i (T “V “ v 



and 



x 2 (t) = x 2 + a 2 y 2^ T-T v^’ 



a - ) b ,y 1 



y 1 ( T ) = exp{b 1 x 2 (T-T v ) + 2 -- 1 — 2 - (T-T v ) 2 }, 

/ X f 

y 2 (T) =y 2> 

dpi 



" a U y l P l 



with P 1 ( T=T V > = P x » 



dt 

dp^ v 

— = - bl y iqi with P 2 (t=T v ) = pj, 

dq l X 1 v 

dT = - a n -f - b ! x 2 q i with q i (T= V = V 



^2 v 

?r = ' a 2 p 2 With q 2 (T=T v ) = q 2- 



We have not been able to develop solutions in terms of "elementary" functions 
to the equations for x^, p^, p^, q^ * 
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Table I. (cont.) - 4 



5. P 



II 
A2 * 



( U * (T) = 1 n 

< for T £ T £ T 

‘v*(T)=l SL 



x is the smallest positive root of S (t=t ) = 0, where 



c ( , C SLII . , SLII 

S V (T) = S v + a 2 b l y l 



X ll( a il f SLII f Y l 

fe 1 k 



‘lfl , II. a il b l X 2 (T_T SL ) 

— f| (t_t sl } " 7. ' f ' 6 



b l X 2 



(b l X 2 )2 



}• 



Again, an upper bound on is given by a i2^ a 2^1^ * ^ as ^ een s ^ own 

that S (x) > 0 for x** < t £ t . Also, on P** we have 



‘SL 



A2 



SLII , f, II, 
b i x 2 (t-T SL ) 



, \ SLII / SLII, II. , a li y i f b l X 2 (T_T SL ) ~1\ 
X 1 (T) = X 1 exp i a 12 y 2 (T_T SL ) + — f — L 6 • 1 J/> 



b l X 2 



x 2 (t) = x 2> 

, f , II, 

, . SLII b l X 2 (X T SL )s 

y i( T ) m y 1 e 



, v SLII 
y 2 ( x ) = y 2 



and 



SLII , f, II, 
b 1 X T ^ T_T SL 



, , SLII ( SLII , II a n y l r b lV T X SL ; "]l 
Pl (T) = P 1 exp (- a l2 y 2 (T-T SL ) T f [f ~Ui 



b l x 2 



with p 



SLII 



1 f X l| 



SLII I f 
X 1 y l 



, , SLII , SLII 
P 2 (t) = p 2 - 



f 




f 




X 1)J a ll + 1 


r slii 


f y ll 


a 

n] 


yf (b.x^) 2 


L q i 


f 

*1 


b x f J 



3]- 



n, 

t sl ) 



, f / n, 

a ll b l x 2 (T_T SL ) 
e 



<V2> 2 



}• 



, f, II, 
, , SLII - b l X 2 (T - T SL ) 
q i( x ) = q x e 



11 



{ 1 -e _b l^ (T " T ^ ) }, 



b„ x 



12 



( v SLII 
q 2 (T) = q 2 " a 12 



f 

V 



, II, 
^sl)' 
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x / u*(t) = 1 

> . P. 0 : \ for t & x & t 

A2 ' V*(T) « 0 



Table I. (cont . ) - 5 
I 



Results are similar to those for P A , above in 4. 

A1 



7. P 






u*(x) = 1 
v*(x) = 0 



* 1 

for T i T T 
V u 



Results are similar to those for P A , above in 4. 

A1 

J / u*(t) = b 2 /(b 1 +b 2 ) • (l-q 2 y 2 /(p 2 x 2 >) = u* 



'• P S2 ; { 



v*(x) = 0 



hold at t = t 



Also, on P ^2 we have^ 

_ 

dT " a il X i y i 
dx 2 

d — = a 2 y 2 

dy i * 

d r - u s b i x 2 y i 

d V 2 * 

— = (l-u s )b 2 x 2 y 2 

and 



d Pi 

dT = " a n y i p i 
dp 2 

dT” = _b i y i q i 



dqj 

dx 

^2 

dx 



"ll f S u l"2 H l 



X 2 V 2 



S 2 2 4 2 



: x 2 )) = 


* 

u s 


^ “I 




for 


X i X ^ X 

v u 


for 


T > X . 
V 


In order for a 


•k 


■- x the following condition must 

V 


* 


• * 




- q 2 (T v 


» y 2 ( V 






* 


V* 


with 


x 1 (t =t v ) 


= x i • 




* 


V* 


with 


x 2 (t=T v ) 


- X 2 > 




* 


V* 


wi th 




= y l > 




* 


V* 


with 


y 2 (x=T v ) 


= y 2 . 




* 


V* 


with 


P l (T=T V ) 


= p l ’ 




* 


V* 


with 


p 2 (T~T v ) 


= P 2 ’ 




* 


V* 


with 


q i (T=T v ) 


= q l • 




* 


V* 


with 


q 2 ( T=T v ) 


= q 2 . 



+ A further discussion of the continuity of the adjoint variables is to be found 
in Section 5.2.4. below. 
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Table I. (cont.) - 6 



9. P 



II 

B1 



■ { 



u*(T) =1 y 

for 0 ^ t ^ t with —p < 
v*(x) = 1 U y 2 



y 2 



is the smallest positive root of 



K a H y i 
b i ' —r~ 



- a 12 b 2Vu + f 



a ll y f b lVu 



X, 



= 0. 



3t 



It should be noted that -r— ^ - > 0, where r = yf/y». It may be shown that for 

or 12 



0 < 



T<(^) 



a l2 y 2 



T < T < T > 
_f u u 



where the determination of t is given in 1. We also have that 

u 

T^(r^) < T u ( r 2 ) f° r r i < r 2 ^ x i anc * x 2 ^eld constant). The solutions 
to the state and adjoint equations are the same as those for PgJ given above, 



Let S (t=t ) = S U , p (t=t ) = pi 1 , 
v u v 1 u r l 



etc. 
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Table I. (cont.) - 7 



IX , u*(t) = 0 

10 • p r9 : \ for x Si x ^ x 

B2 ' v*(x) = 1 u v 



It follows that for all t > x we have S (t) < o and — (t) < 

U U V 2 

a l2 b 2' ^ a ll B l^ ' T v t * ie sma H est positive root of S v (t=x ) = 0, where 
S (t) is given by 



S v (T) = S v + a 2 b 2 y 2 



(X, 



y f 



(b/j' 



- 


f 








u 


f y i] 


+ 312 


(T-T u ) - 


a l2 


q 2 


l x f J 

X 1 


b 2 x 2- 


< b 2 x 2 )2 



exp[b 2 X2(x-x u )]|. 



Also, on we have 



f . f, 

n .. a n 9 y 9 b x (t-x ) 

x l (,) = x i ex p^ a 11 y 1 ( T-T u > + - t« 2 2 -1]>, 

b 2 x 2 

x 2 (x) = x 2 , 



y l (T) = y i» 



y 2 (T) = y 2 ex P^b> 2 x 2 ( T-T u ) ^ * 



and 



f , f, . 

u u b^x^(x T ) 

Pl( T > “ P-j^ exp{-a 11 y 1 (T-T u ) p [e 2 

b 2 x 2 



■1] } with p“ = 





f 


- /^X_ >1 


1 


1 


u 


f 


X 1 


V 



/ \ u , f 
P 2 M - P 2 - b 2 y 2 



/•x. 






12 



( b 2 X 2^ 2 



r> 


f 


— 


u 


r y i^ 


+ ^2 


q 2 


f 


+ , f 




^x. 


b _x_ 




1 


2 2-* 



(x-x u ) - 



12 



b„x (x~x ) 
2 2 u 



(b 2 x 2 )2 



qi(T) - q“ - a u 



f 

<X^ 



(X 



q 2 CO = 



7, 



12 



b^x. 



u 




_ 


f 


a i 0 

12 1 , u 


f y il 


U f q 2 


L f . 


b „x - 


l x. 


2 2 


1 



i - b 2 x 2 <T - T u>> 
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Table I. (cont.) - 8 



11 . 



II / u *(t> 0 II 

P : { for * x & x 

63 < v*(i) - 1 SL 



t is the smallest positive root of S (t=t ) = 0, where 
v v v 



e , v C SLII , . SLIlf X l 

S v (t) = S v + a 2 b 2 y 2 -r 




12 



y l ' ^ b 2 x 2 ^ 



*12 






, II , 
(t " i sl ) - 



p 


f 


SLII 


f y l] 


q 2 


L f J 

x i 




a 12 


t 



f , II, 
2 x 2 < T - T SL ) 



( b 2 x 2 ) 2 



Again, an upper bound on x r is given by a (a 2 b^) . It may be shown that 



v 



II 



II 



S (x) < 0 for all x > x CT . Also, on P we have 



SL 



B3 



, , SLII , SLII , II, . a 12 y 2 
x 1 (x) = x 1 exp{a 11 y 1 ( t-t S l^ + T 

b 2 X 2 

x 2 (x) = x 2 . 



SLII . £, II, 

b 2 x 2 ( x - x C T ) 



[e 



SL' 



- 1 ]}, 



, , SLII 
y i ("0 = y ± 



y 2 ( T ) = y ^ 11 exp{b 2 x 2 (x-Xg^) }, 



and 



SLII f, II. 

, , SLII f SLII , II, a 12 y 2 r b 2 X 2 (T T SL' , ,, . , SLII 

p l ( T ) = P x exp{-a 11 y 1 ( T-T SIj ) f — t e — 1 3 > with p 1 



. , SLII . SLII 
P 2 ( t) = P 2 - b 2 y 2 



f 




b 2 X 2 


f 








f x ii 


J ^2 + 


SLII 


r y ii 


+ a 12l 


(t-t sl 3 


a 12 


f 

v 


(b 2 x 2 )2 


q 2 


l x f J 

X 1 


b 2 X 2- 


(b 2 x 2 ) 2 



1 


f*ll 


SLII 


f 


X 1 


N 

1— 1 



b 2 x 2 (T ‘ T SL ) 



}• 



, , SLII 
“l < ' l) = q x - a u 



rX 



7 i 



1 , II, 
?P SL 3 ’ 



< 2 « - ft 



{- 



12 



b 2 x 2 



f 1 



12 , SLII 
~ +q 2 



(V, 



>2 X 2 



f, 

V - 



- b 2 X 2 (T ‘ T SL ) 
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Table I. (cont.) - 9 



12. P 



B4 



•{ 



u*(x) = 0 
V*(T) = 0 



for t it 
v 



It may be shown that S U ( T ) < 0 and S^CO < 0 for all t > t 



on P . we have 

DH 



V* r V* * i 

x i (T) = x i expta n y i 



4 r > 2 - 



2a 



2 v* 



2 , v *. 2 



v* 



b 2 y 2 coth ( -A ( T-T v ) +B ) for “ 2 “ ( x 2 > > a 2 y 2 » 



x 2 ( T ) = ^ xf/U x 2 *( t - t *)) 



t 2 , v*. 2 v* 

for -y ( x 2 > = a 2 y 2 ’ 



7 ^ 



2 V* v* 9 "k ? V* 9 V* 

b - y 2 - ( x 2 ) 2 tan(C(T-T v )+D) for y (x 2 ) 2 < a 0 y 0 , 



2 y 2 



V* 

y 1 (T) = y l ’ 



/b (x 2 ) 2 v *\ * b 2 . v*. 2 v* 

y 2a y 2 J/sinh 2 (-A(T-T v )+B) for — <x 2 ) 2 > a 2 y 2 , 



y 2 (t) = \ y T /{l - T X 2 (T - T v )} ' 



, 2 , v*. 2 v* 

for — (x 2 ) 2 = a 2 y 2 , 



(y 



V* 



^2 v* , o * ^2 v* 9 v* 

2 ^” ( x 2 ) )/cos 2 (C(t-T v )+D) for — ( x 2 ) 2 < , 



. . b 2 A v* 2 2a 2 v* 

where A = — / (x 2 )* - — — y 2 , 



Also , 
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Table I. (cont.) - 10 



12. P , : (concluded) 

d4 



and 



v* 



B = coth 



-11 



r f v* 2 2a 2 v*y 

U 2 ) - TT y 2 



c = 



/f!2 v* _ v* 2 

2 v b y 2 ^ 2 ; 



v* 



D = tan 



-11 






[/S 

V b_ 



^ a 2 v* , v* 2 

b y 2 ( x 2 ^ 



p x oo 



V* r V*. * , . , V* 

P 1 expi-a^j^ (x-x^)} with p x 






V* 



X 1 y l 



d -h- k 

dx b 2 y 2 q 2 



. , / *\ V * 
with p 2 (x=x v ) = p 2 , 



q x (t) 



v* 

= q. - a 



11 



f 



V 



(t-T*), 



^2 

dx a 2 P 2 " b 2 X 2 q 2 



* V* 

with q 2 (T=T v ) = q 2 . 



We have not been able to develop solutions in terms of Elementary' 1 
to the equations for p 2 and q 2 * 



functions 
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Table I. (concluded) - 11 



13. Pg 2 and P 



1 , u*(x) = 0 

B3 " { for t ^ x 



v*(t) = 0 



Results are similar to those for P„, above in 12. 

B4 



u*(t) = 0 



14 



' P B5= U.o X SL 



^ T 



Results are similar to those for P_, above in 12. 

B4 



15 



I / u *(t) 1 i i 

• P./ '• < for T ^ T ^ T 

A4 1 j. / n SL u 

v v*(t) = 0 



Results are similar to those for P A1 above in 4, 

A1 



5.2.3. Numerical Examples . 
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In this section the results of some numerical computations are presented 

and discussed. A computer program which calculates the information given in 

Table I was developed in FORTRAN for the IBM 360 computer. I would like to 

thank CPT Jeffrey L. Ellis, U.S. Army for his efforts on this project. This 

computer program outputs a plot representing the field of extremals (see 

Figures 5 and 6 below). The closed-form analytic results presented in Table I 

are used whenever possible. Approximate numerical solutions to transcendental 

*)(• 

equations (for the determination of, for example, t^, t , etc.) are developed 
by the well-known Newton-Raphson method. In those cases for which closed-form 
solutions are not available to the state and adjoint equations, a standard 
fourth order Runge-Kutta numerical integration method is used. A time step, 

At, was used in these numerical integrations which yielded agreement to the 
fifth place to the right of the decimal place in test cases in which the approx- 
imate numerical solution could be compared with the exact solution. 

Parameter sets for the numerical examples considered in this report are 
given in Table II. In the problem (64) at hand we may consider time to be an 
additional state variable so that the state space in five dimensional, i.e. 
the state variables are t, x^, x^, y^, and y^. Thus, unfortunately, we 
cannot graphically depict the field of extremal trajectories, since this would 
involve curves in a five-dimensional space. 



Table II. 

Parameter 

Set a ll 

1 0.003 



Parameter Sets Used to Generate Numerical 
Results Shown in Figures 5 and 6. 




0.006 0.01 0.004 0.005 4.0 8.0 8.964 



2 



0.003 



0.006 0.01 0.004 0.005 4.0 8.0 11.597 
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The most illuminating plot for gaining insight into the structure of the 
optimal fire support strategies for (64) is that of extremal trajectories in 
terms of y /y£ versus backwards time, x. This is shown for parameter set 
1 in Figure 5. The corresponding strategic variable values for X and Y 
(i.e. u* and v*) along each extremal are also given. Other plots have 
been considered, but they provide little, if any, additional insight. 

The most significant features of the field of extremals shown in Figure 
5 are the two U-singular "surfaces": there is one in £,£ - £,cj^ space in 

V-phase I and one in £-space in V-phase II. In each phase, X uses the 
strategy U = 1 above the singular "surface" and the strategy U = below it. 
Similar to our discussion in [27], the singular surfaces are present in the 
field of optimal trajectories so that the X artillery avoids "overkilling" 
either or Y^. This is obvious when one, for example, considers 




Thus, the rate of destruction of Y^ per unit of X artillery decreases over 

time as the Y^ force level decreases (see [26] and [27]). 

Results for parameter set 2 are shown in Figure 6. There is a void 

(see p. 169 and also p. 187 of [15]) in the field of extremals. This is because 

* 

in backwards time at the end x^ of the U~singular subarc in V-phase II, we 

& 4 - 

would have u c (x ) (as given in Figure 3) equal to 1.054 if the adjoint 
o v 

£ 

variables were continuous at x • The following theorem further explains this 



situation. 
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a 12 b 2 

a ll b l 






Backwards Time, t (minutes) 



(t=T) 



Figure 5. Plot of — versus Backwards Time, t, for 

y 2 

Field of Extremals for Parameter Set 1. 




x , for Parameter Set 2. 



versus Backwards Time, 
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THEOREM 6: There can be no U-singular subarc beginning in 

backwards time at x + with a^b^, = a 10 b_y n for 

v 11 rl 12 2 2 

b l P 2 (T v )x 2 + b 2 q 2 (r v )j, 2 ' °' 

When a U-*singular subarc begins at x* with 

a ll^l^l = a l2^2 y 2 > t ^ iere no discontinuity in the 
adjoint variables at t = (i.e. a = 0 in (99)). 



PROOF: Immediate by (89) and (99). 



Q.E.D. 



Additionally, Theorem 7 gives the extremal transitions in X f s strategy possible 

* 

from the U-*singular surface in V-phase II as we work backwards from x , Thus, 

since b n p 0 (x )x n < b 0 (-q 0 (x ))y n for parameter set 2, a void would exist in 
12 v 2 2 2 v 2 

k 

the field of extremals if the adjoint variables were continuous at x . 

J v 



THEOREM 7: Assume 


that there is 


no discontinuity in the adjoint 


variables at x 


= X 

V 


with a 11 b 1 y 1 = a i2 b 2 y 2* 1111611 


I. if 


b l P 2 (T v )x 


2 " b 


2^ _q 2^ T v^ y 2’ then We Can ° nly 


have 


u*(x) =0 


for 


t £ (t^jT^+ 6) where 6 > 0, 


II. if 


b iP 2 < T v >* 


2 * b 


2 ( -< ?2 ( T v ) ^ y 2 ’ then We Can have 




1 


' (a) 


0, 




U*(T) = \ 




(l- q 2 y 2 ^P2 X 2^ * b 2 /( b 1 +b 2^ 5 




1 


v (c) 


1, 


for x 


C: (x , x +6) where 6 > 0. 

V V 



PROOF: (a) When we are on the singular surface in V-phase II at x = x 



*- 



then by (84) and (85) and the continuity of the dual variables we have 



S (x=x + ) = S (x=x + ) = 0, 
u v u v 



and 
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s u ( t=t+ ) = a 2 b 2 (b 1 +b 2 )p 2 x 2 y 2 {u*(T+)- 



CN 


fl V2l 


lVb 2 ; 


l p 2 x 2 J 



} , 



(109) 



o dS 

where S denotes 

u dx 



u 



(b) Considering a Taylor series expansion about x = x , we have by 

v 



+ *+ 

the above for x £ x = x 

v v 



S u (T) - 



* O 

(T-T v ) 2 oo 



S u (T), 



( 110 ) 



where x € (t^,x). 



(c) When u*(x) = 0 for x £ (x*,x*+S), then 

v v 



S u (T ‘ T v +) * ■ a 2 (b 2 )2p 2 x 2 y 2 



1 - 



q 2 y 2^ 



P 2 V 



< 0, 



so that 3 5 > 0 such that S (x) < 0 for all x £ (x x +6 ) . Thus, we 

1 u v’ v 1 

can always have u* = 0 as we work backwards in V-phase I from the U-singular 
subarc in V-phase II. 



(d) Now let b 1 P 2 ^ T v ^ x 2 ^ b 2^~ q 2^ T v^ Y 2’ By ( 88 ^ > the U-sin 8 ular 
control in V-phase I u = (l-q 2 y 2 / (P 2 X 2 ) ) •b 2 /(b 1 +b 2 ) £ 1. Thus, the 

*+ ° ° * 4 - 

U-singular subarc is possible. When U *( T V ) = 1> then S^(x=x^ ) ^ 0 by 
(109). When inequality holds, it follows that 3 ^ > 0 such that S^(x) > 0 

•k Jc * 

for all x £ (x ,x 4-6-). Clearly, we cannot have u* = 1 if b-p 0 (x )x 0 < 
v v 1 12 v 2 

b 2 (“ q 2 ( T v^ y 2‘ Q.E.D. 

The same analysis as used in the proof of Theorem 7 applies on a 
U-sintular subarc in V-phase I when v* = 0 . As long as (89) holds, one 
has three options similar to those of part II of Theorem 7. 
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5.2.4. Filling in a Void . 

We have emphasized that H, £(t) , and cj^(t) are continuous function 
of time except at manifolds of discontinuity of both U* and V* (see Section 
4.3 above). From Theorem 7 it follows that a void must exist in the field of 

k : k 

extremals when these functions are continuous and ) x £ < ^2^”^2 ^ T v^^2’ 

k 

At moreover, v* changes (as we progress backwards in time) from 1 to 0 

and u* from b^/ (b^+b^) to a different value. Thus, we have a manifold of 
discontinuity of both U* and V*. Moreover, considering results given above, 
it is readily shown that u*(x) remains for increasing T (i.e. backwards 
time) equal to zero once it changes to zero. Then from Theorems 6 and 7 it 

k _ ^ 

follows that for bj^^v ) x 2 < ^^2 t ^ e var i a bles must be 

* I* 

discontinuous to fill in the void, and we must have u*(t) - 1 for T < t < T 

v u 

Furthermore, considering Figure 6 and considerations n in the large," the mani- 
fold of discontinuity must lie on the V-transition surface. 

Thus, we have established that for a ^b^y^ = a ^2^2 y 2 we b ave 



/ U *( T V ) = b 2 /(b 1 +b 2 ) , 
' v*(t ) =1 



and 



( u *< T v > = 
v*( 7+ ) = 0. 



v 



3t 



( 111 ) 



3t 



and 



It remains to determine the function x (x,y) of (91) so that — — ^ - 

v ~ 8x 3^ 

may be computed, and the jumps in H, £, and subsequently determined (see 

(92) through (95)). It should be clear that it is impossible to explicitly 
determine t^Oc,^). However, by computation of five points on the V-transition 
surface, the desired partial derivatives may be estimated by using linear 



approximations to the appropriate directional derivatives and solving a system 
of four linear equations in four unknowns. For parameter set 2 (as the reference 
case), this yielded the following estimates 
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3t 


3x 


- 0.0000, 


V 


8x l 


3x 2 


3t 


8t 



3 yi 



= -0.0167, 



3 y^ 



= -0.295, 



= -0.0331. 



It is, therefore, convenient to re-write the jump conditions across 
the manifold of discontinuity of both U* and V*. 

3x 



p i <T v> ■ p i ( V- 



MO ■ - p 



+ 3t 
W ’ <1 l (T v ) - p 3^ - oa ll b l- 



(113) 



+ _ 3t 

q 2 ( V 'VV - P 8^ +0a 12 b 2’ 

where p and a are related by (98). In this case the jumps (99) and (100) 
in the switching functions simplify to 



S U (T V ) * °kl <b l )2y l + a 12 (b 2 )2y 2 + 



a ll <b l )2y l x 2 



3x 



3x 



|Vl 3y_ 



- b„y 



r 2 3y 2 ) 



3x 



l-a 0 y 



2 J 2 3x 



3x 

v , v 

b, y n x -r — 

112 3y p 



(114) 



and 



3x 

a ll a 2 (b l )2y l X 2 3x7 ° 

S (x + ) = 

V V 



8x 



l-a 0 y 



2 2 3x„ 



9X vl 

- b l y l x 2 3^] 



(115) 



^ I 

Since v*(x ) = 0, we must have S^(x v ) ^ 0 so that (112) and (115) 
yield that o ^ 0. It should be clear that a = 0 if and only if H, £, and 



* 

q are continuous at x . 
^ v 



that we must have 



*+ 

For a > 0, the condition that U *( T V ) = 1 yields 



*+ 

S (T ) 

U V 



> o, 



(116) 



where ^ U ( T V ) i s gl ven by (114). Although it cannot in general be guaranteed 
that (116) will always hold when a void in the field of extremals such as that 
shown in Figure 6 exists, it should be clear that it must if the problem (64) 
is to have a solution. The author conjectures that this is true. It is readily 
shown that when (116) holds, we have 



S 

u 




) > 0 , 



S 



u 




) < 0 , 



and 



S 

v 




) < 0 . 



(117) 



As we have shown in our past work (see Appendix A) , the appropriate 

value for a is determined by "considerations in the large": the structure 

of the entire field of extremals determines the value of this parameter. In 

I* 

Figure 7, we let x^ denote the backwards time at which the U-singular subarc 

I* 

is entered in V-Phase I. Corresponding to x^ is a*, which yields the 
first and second conditions (86) and (87) (with u* ^ 1) for a U-singular 

I* ^ 

subarc with V* = 0 at x > x . For 0 < o < a*, one uses u*(x) = 1 for 

u v 

*+ i X . 

x < x < x and then u*(x) = 0 for x > x . For a > a* , the U-switching 
v u u 

* 

function S^(x) never changes sign so that u*(x) = 1 for all x > x^. Thus, 
by manipulation of a, one may fill in the void in the field of extremals in 
V-Phase I. The resulting field of extremals is shown in Figure 7. 



5.2.4. The Case of Negligible Small Arms Effectiveness . 

It seems appropriate to consider what happens to the solution to the 
problem at hand as the (relative) effectiveness of (small arms) fire 

becomes negligible, i.e. as a^ 0. Let us consider (either) Figure 5 (or 
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Figure 



7. 



Filled-In Void in Field of Extremals 
For Parameter Set 2. 
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Figure 7). The O-singular '’surface' 1 in V-Phase II has equation = 

a l2^2^ a ll^l^ * Thus, as a^^ ■> 0 with the other parameters being held con- 
stant, this singular "surface" appears higher and higher on the y^y^ axas 
in Figure 5. In the limit, the singular surface does not appear in the finite 
part of the plane. Thus, we have shown that an optimal strategy in which a 
side divides the fire of its supporting weapon system between the enemy's 
primary (infantry) and supporting systems can only occur when the enemy’s 
infantry has some effectiveness (in the sense of a non-zero Lanchester attri- 
tion-rate coefficient) against his infantry. 

6 . Discussion . 

In this section we discuss what we have learned about the dependence 
of the structure of optimal fire-support strategies upon model form from our 
examination of the problems at hand. We studied two aspects in this appendix: 
(1) the relationship between the criterion functional and the structure of 
optimal fire-support strategies in Kawara’s fire support differential game 
[16], and (2) a different attrition form for casualties produced by the 
attacker’s fire support. 

We first examined for what class of criterion functionals (terminal 
payoff only) the optimal fire-support strategies (at least for the case in which 
neither side’s supporting units are all lost) are independent of force levels. 
This was done for the same combat dynamics as considered by Kawara [16]. In 
this examination, we introduced a general approach for such determinations. 

We determined the most general form of terminal payoff with only infantries 
considered in the criterion functional for which optimal fire-support strate- 
gies are independent (as qualified above) of force levels. The criterion 
functional considered by Kawara [16] is a special case of this general form 



D-67 



of terminal payoff. Thus, Kawara [16] essentially used the only objective 
function that yields the results reported by him. For any other form of 
terminal payoff, the (forwards) time at which a side changes from concentra- 
tion of all fire on the enemy’s supporting units (counter-battery fire) to 
his primary ones in following an optimal strategy will depend upon the force 
levels. Finally, we showed that there is no ’’simple" type of terminal payoff 
which includes both sides’ infantry and artillery and yields optimal fire- 
support strategies which are independent of force levels. 

We finally examined the dependence of optimal fire-support strategies 
upon the form of the combat attrition model by considering slightly different 
combat dynamics than those considered by Kawara [16]: we let the attacker’s 

(X’s) artillery produce "linear- law" attrition against both the defender’s 
artillery as well as infantry. As contrasted with the optimal strategies in 
Kawara ’s problem [16] of always concentrating all artillery fire on first enemy 
artillery and then enemy infantry, for our problem (64) the optimal strategy 
of one combatant (X) depends directly on the enemy’s force levels and is no 
longer to always concentrate all fire on either the enemy’s primary or secondary 
weapon system. The solution to (64) is characterized by the presence of singu- 
lar surfaces (a different one for each V-Phase of battle) . When the system 
state reaches one of these surfaces, X follows a strategy of dividing his 
artillery fire between enemy infantry and artillery in order to avoid "overkill." 
Another aspect of optimal fire-support strategies not present in Kawara’s 
problem [16] which occurs in the problem at hand (64) is that X’s strategy 
depends upon the optimal fire-support strategy of Y. In fact, X must react 
instantaneously to changes in Y’s optimal strategy. 
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It should be pointed out that the development of even a partial 
solution to (64) involved solution phenomena not previously encountered in a 
Lanches ter- type differential game: the adjoint (or dual) variables^ are 

discontinuous on a manifold of discontinuity of both U* and V*. This 
manifold of discontinuity exists for a certain range of parameter values in 
the solution to the problem at hand. Furthermore, there is a military inter- 
pretation to this manifold of discontinuity: if concentrates fire on 

and on Y^, then when Y^ changes to concentrating all fire on 

X^, X must re-evaluate the worth of a Y 2 unit because it now has a direct 
influence on the payoff. Such discontinuities in the dual variables are unique 
to differential games (see Appendix A). 

It should also be pointed out that the presence of the singular surfaces 
in the solution to (64) is independent of the form of the criterion functional 
(here terminal payoff) and depends only on the combat dynamics. For purposes 
of comparison, we considered the same payoff as considered by Kawara [16] . 

We also showed that the singular surfaces can only be present in the solution 
when the defender's infantry (Y^) has a (non-zero) casualty producing capa- 
bility against X^ . 

The problem (64) considered in this appendix has certain similarities 
to the "War of Attrition and Attack: Second Version" considered by R. Issacs 

(see pp. 330-335 of [15]). We have, however, developed a much more complete 
solution to the problem at hand. Although this problem (64) possesses some 
similarities to the Lanchester-type optimal control problem studied by us in 
[26], its solution has turned out to be much more complex. It should be pointed 



t 



The reader should recall that these represent the marginal values of force types, 



3V 

i ' e - P 2 (t) * 3^00 ” here V 
game (see [10], [15]). 



V(t,£,^) denotes the value of the differential 
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out, however, that our developments in this appendix have been materially 
aided by intuition developed through the study of the simpler, one-sided 
problem (see [27] for a further discussion). 

As a result of our investigation here, we hope that a better under- 
standing of optimal fire-support strategies has been developed. Our work here 
shows that the nature of the functional form of the artillery casualty rates 
essentially determines the most significant aspects of the structure of 
optimal strategies. Thus, examination of these optimization problems shows 
the importance of determining the appropriate (Lanchester-type) model of 
combat dynamics. 

Finally, it should be pointed out that computational work on the problem 
(64) has not been completed at the time of this report. Due to the theoretical 
complexities (especially the discontinuity in the dual variables), student 
project work on numerical examples was not completed. We would propose to 
ONR the completion of this computational work as a future research task. 
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